THE RELATIVISTIC SCOTT CORRECTION FOR ATOMS AND 

MOLECULES 
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Abstract. We prove the first correction to the leading Thomas- Fermi energy for the 
ground state energy of atoms and molecules in a model where the kinetic energy of the 
electrons is treated relativistically. The leading Thomas- Fermi energy, established in [25], 
as well as the correction given here are of semi-classical nature. Our result on atoms and 
molecules is proved from a general semi-classical estimate for relativistic operators with 
potentials with Coulomb-like singularities. This semi-classical estimate is obtained using 
the coherent state calculus introduced in [36]. The paper contains a unified treatment of 
the relativistic as well as the non-relativistic case. 
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1. Introduction and main results 

Our goal in this paper is to study how relativistic effects influence the energies of atoms 
and molecules. More specifically, we are aiming at proving a relativistic analog of the cele- 
brated Scott correction [29, 16, 13, 15, 30, 31, 32, 36]. At present there is no mathematically 
well-defined fully relativistically invariant theory of atoms and molecules. We will here con- 
sider a simplified model, which shows the relevant qualitative features of relativistic effects. 
In this model, these effects are introduced by treating the kinetic energy of electrons of mass 
m by the operator V— ^^c^A -|- m^c^ — mc^ instead of the standard non-relativistic Laplace 
operator — /i^A/2m. Here c refers to the speed of light and h is Planck's constant. It is the 
simplest of a class of models that attempts to include relativistic effects; see [12, 23]. Al- 
though this model does not give accurate numerical agreement with observations it is from 
a qualitative point of view quite realistic. 
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One of the qualitative features that our model shares with all other relativistic models is 
the instability of large atoms or molecules. The natural parameter to measure relativistic 
corrections in atoms and molecules is the dimcnsionlcss fine-structure constant a = /he, 
where e is the electron charge. As we will explain below, if Za is too large {Z is the atomic 
number) then atoms are unstable. In our model the critical value of Za is too small compared 
with experimental results and with what is assumed to be the correct critical value, namely 



Our main interest here is the behavior of the total ground state energy of large atoms and 
molecules. Because of the relativistic instability problem mentioned above one cannot simply 
consider the limit of large atomic number Z. One is forced to look at the simultaneous limit 
of small fine-structure constant a in such a way that the product Za remains bounded. Of 
course, a has a fixed value which experimentally is approximately 1/137. Thus considering 
the limit a tending to zero is strictly speaking not physically correct. Likewise, considering 
the limit of Z tending to infinity is in contradiction with the fact that the experimentally 
observed values of Z are bounded (by 92 for the stable atoms). Studying the limit Z 
GO and a — >^ with Za bounded allows us to make a precise mathematical statement 
about the asymptotics. There is numerical evidence that the asymptotics is indeed a good 
approximation to the total ground state energy for the physical values of Z and a. 

The first to, at least heuristically, suggest to consider Za as a separate parameter in the 
limit Z — > oo was Schwinger [27]. In this original paper, Schwinger finds discrepancies of his 
estimates of relativistic corrections with numerical evidence. Later [8], more corrections are 
taken into account and excellent agreement is found. This accuracy however goes beyond a 
rigorous mathematical treatment. We content ourselves with giving a rigorous treatment of 
the simplified model with the correct qualitative behavior. 

The first rigorous treatment of the limit Z — > oo with Za bounded was given by one of 
us is the paper [25], where the leading asymptotics of the ground state energy was found. It 
turns out it does not depend on Za. The goal of the present paper is the first correction to 
the leading asymptotics, i.e., the Scott correction and, in particular, to show that it depends 
on Za. The work in [25] was generalized to another relativistic model in [4]. 

We now introduce the molecular many-body Hamiltonian we consider in this paper. Let 
e and m denote the electric charge and mass of an electron. Let Ze = (Zie, . . . , Zm©), 
where Zi, . . . , Zm > 0, be the charges of the M nuclei. We consider the Born-Oppenheimer 
formulation where these nuclei are at fixed positions R = . . . ,Rm) £ M"^^^. We have 
N electrons. As explained above the relativistic kinetic energy of the j-th electron is equal 
to y^—h'^c'^Aj + m^c^ — mc^, where Aj is the Laplacian with respect to the j-th electron 
coordinate yj & M.^, j = 1, . . . ,N. The potential energy of the electrons is composed of the 
attraction to the nuclei, 



Za = 1. 




(1) 



and the electron-electron repulsion, 




The total energy of the electrons is described by the Hamiltonian, 




Vi - Vj 
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Let us now introduce the fundamental constants. Namely, let a = /i^/me^ be the Bohr radius, 
and Roo = ^m.e'^/h^ Rydberg's constant. Then by a change of coordinates yj — Xj = Uj/a, 
we see that 

(2i?oo)"^^^rei =: H{Z, R; a) = H{Zi, . . . , Zm, Ri, ■ ■ ■ , Rm\ ot) 



L • -II JUt JLi 1 I 

j=l \<i<j<N ' * 

where again a is the fine-structure constant. For a = the kinetic energy of the j-th electron 

Here we have set N = Z = X]yt=i so that the molecule is neutral. In particular, this 
means that Z must be an integer. Prom now on we study the operator H{Z,'R;a). This 
operator acts as an unbounded operator on the anti-symmetric tensor product, /\^ L^(M^ X 
{— 1, 1}), where ±1 refers to the spin variables. We arc interested in the ground state energy 

E{Z, R; a) = mfa{H{Z, R; a)) , 

and, in particular, in an asymptotic expansion of this when Z oo. 

The ground state energy E{Z, R; a) is finite if max^j^fca} < 2/-K, but E(Z, R; a) = — oo 
if maxjtjZjfcQ;} > 2/7r (see [12, 23])-'^. This is the relativistic instability discussed above. 
Therefore we must require the atomic numbers to be smaller than or equal to 2/(7ra) which 
is approximately 87. This is of course in contradiction with the experimental fact that larger 
stable atoms exist and is one reason why our model can only be qualitatively correct. (We 
want to emphasise that the instability we are discussing here is not the nuclear instability 
causing atoms larger than atomic number 92 to be unstable. The relativistic instability we 
discuss here is only believed to manifest itself for atomic numbers greater than 137.) 

The true energy of the molecule should include the nuclear-nuclear repulsion. Since the 
nuclei arc considered fixed here the nuclear-nuclear repulsion is simply a constant which we 
have omitted. 

As discussed above the leading asymptotics of E{Z, R; a) will be independent of the rel- 
ativistic parameter a. It will be given by what is called Thomas-Fermi theory. The seminal 
contribution by Lieb and Simon [21] was to put Thomas-Fermi theory on a solid mathemat- 
ical foundation and to prove that in the non-relativistic case the Thomas-Fermi energy of a 
molecule is indeed the correct leading asymptotic energy for the true ground state energy 
as Z — ^ oo. This is the result that was generalized to our relativistic model in [25]. 

The main result of this paper is the following asymptotic result on the ground state 
energy. 

Theorem 1 (Relativistic Scott correction). Let z = {zi, . . . , zm) with zi,...,zm > 0, 
"^k^iZk = I, and r = (ri,...,rM) G with miufe^^lrfe — r^l > ro for some ro > 

be given. Define Z = {Zi, . . . , Zm) = Zz and R = Z~^/^r. Then there exist a constant 
E'^^{z,r) and a universal (independent ofz, r and M) continuous, non-increasing function 
S : [0, 2/7r] M with 5(0) = 1/4 such that as Z = Ylik=i'^k ~* oo and a ^ with 
maxfc{Z/;Q!} <2/iT we have 

M 

E{Z, R; a) = Z^^E'^^iz, r) + 2 ^ ^|<S(Zfea) + ©(Z^-VSO) _ (2) 

fe=i 



^Hcrc, and in the sequel, operators are defined as the Friodrichs extension for the corresponding form 
sum, originally defined on C^f-functions (here, for instance, C^{M? x { — 1, 1}))- 
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Here the error term means that < CZ^-Vso^ y^/jg^g ^/^g constant C only depends 

on ro and M. As before, ■\/—a~'^A + — = — 5A when a = 0. 

Remark 2. A less detailed version of our result was announced in [35]. 

Remark 3. Several features of our result and its proof should be stressed: 

(i) The constant £^^^(z,r) is determined in Thomas- Fermi theory. 

(ii) The fact that R = Z~^/^r is the relevant scaling of the nuclear coordinates may, as 
we shall see, be understood from Thomas-Fermi theory. 

(iii) A characterization of the function S is given explicitly in Corollary 6 below (see also 
Lemma 25). Its continuity is proved in Theorem 4. 

(iv) The asymptotic result is uniform in the parameters Z^a G [0, 2/7r], k = 1, . . . , M. 

(v) The result contains, as a special case, the non-relativistic situation Z^a = and, in 
particular, the non-relativistic limit is controlled due to the continuity of the function 
S and the uniformity in the parameters Z^a. In order to get the non-relativistic limit 
it is important that all estimates have the correct non-relativistic behavior. This is 
an important issue in this work. Note that in the non-relativistic case the value 
5(0) = 1/4 is explicitly known, whereas this is not the case for any other value. This 
is because the eigenvalues of Hydrogen are explicitly known in the non-relativistic 
case, but not in this relativistic case. The technique to prove a Scott-correction 
without knowing explicitly the eigenvalues for the one-body Hydrogen (like) operator 
was invented by Sobolev [33]. 

(vi) The proof of Theorem 1 does not rely on knowing the non-relativistic case, but treats 
both the relativistic and non-relativistic case simultaneously. 

(vii) The situation near the critical value Z^a = 2/n is understood since the function S 
is continuous up to the critical value 2/7r. This is, however, a less important point 
since we do not know whether the model we study gives a good description near the 
critical value. 

The Scott correction was predicted by Scott [29] to be the first correction to the 
Thomas-Fermi energy. In the non-relativistic setting, this was mathematically established 
by Hughes [13], Siedentop and Weikard [30, 31, 32] for atoms (and by Bach [2] for ions) and 
later by Ivrii and Sigal [15] for molecules. Later a different proof was given by two of us for 
molecules [36]. Based on methods in [15], Balodis Matesanz [3] gave a proof for the Scott 
correction of matter. The Scott correction for operators with magnetic fields was studied by 
Sobolev [33, 34] (in the non-interacting case). 

In [9] , Fefferman and Seco derived rigorously the second correction to Thomas- Fermi the- 
ory for atoms, which is of the order Z^^^. This was predicted by Dirac [7] and Schwinger [28]. 
It is apparently still an open problem to prove this for molecules and to find the relativistic 
correction to this order. 

The main approach to proving the energy asymptotics for large atoms and molecules goes 
back to Lieb and Simon [21] and is to use semi-classical estimates. The Z-scaling makes it 
possible to relate the many-body problem to a one-body spectral problem, which may be 
treated scmi-classically, where the semi-classical parameter is h = Z^^l"^ . Here, several 
techniques have been developed. Lieb and Simon used Dirichlet-Neumann bracketing. This 
is however not refined enough to get beyond the leading term. The Weyl calculus [26] is 
the most advanced and precise method as far as optimal semi-classical error estimates are 
concerned, but it also will not directly give the Scott correction. Ivrii and Sigal [15] used 
Fourier intergral operator techniques to establish the non-relativistic Scott correction for 
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molecules. Hughes and Siedentop-Weikard used methods that were designed particularly 
for spherically symmetric models, i.e., the atomic case. 

A simple method, which is particularly well adapted to many-body problems is that of 
coherent states. It was pioneered by Lieb [16] and Thirring [37] to give very short proofs 
that Thomas-Fermi theory is correct to leading order. It is one of our major contributions 
here to use an improved calculus of coherent states as developed by two of us in [36] to the 
relativistic setting. 

One feature of our work is that we give a general semi-classical estimate for relativistic 
one-body operators for potentials with singularities such as the Thomas-Fermi potential 
(see Theorem 4 below). This is derived by first proving a localised semi-classical estimate 
for potentials with some smoothness (see Theorem 32). The proof here is not much different 
from the one presented in [36] for non-relativistic Schrodinger operators. We do not claim 
that our error estimates arc sharp given the regularity we assume on the potential, but only 
that they are sufficient to prove the Scott correction. In this connection we point out that 
in order to prove the Scott correction it is enough that the error relative to the leading term 
is smaller by more than one power of the semi-classical paramter h. In our case the relative 
error in Theorem 32 is h^^^. 

The relativistic kinetic energy is more cumbersome to work with than the Laplace operator 
and large parts of the rest of our proof from [36] have to be done differently. A main issue 
is to be able to localise into separate regions. Since the relativistic kinetic energy is a non- 
local operator, localisation estimates are more involved than in the non-relativistic setting. 
The philosophy is that localisation errors should behave as if we were working with non- 
relativistic local error terms up to some exponentially small tails (see Theorem 14). 

The proof of the main theorem presented in Section 3 is based on the general semi-classical 
estimate Theorem 4 and the use of a correlation estimate (see Theorem 17) to reduce to the 
one-body problem. 

After we had announced our results in [35], Frank, Siedentop, and Warzel [11] found a 
proof for the atomic case based on the method of Siedentop and Wcikard [30, 31, 32], also 
[10] for the model studied in [4]. This approach seems to be restricted to the spherical case. 
This work does also not, contrary to the present work, make any special treatment of the 
non-relativistic limit or the continuity of the function S. 

1.1. Main semi- classical result. We consider the semi-classical approximation for the 
relativistic operator 



We will consider potentials F : — > M with Coulomb singularities of the form Zfcjx — 
k = 1,... , M, at points ri, . . . , tm € and with charges < zi, . . . , zm < 2/7r. Define 



Tp{-ihV) - V{x) , 



where 





dr{x) = min { |x — r^l | k = 1, . . . , M} , r = (ri, . 
We assume that for some > the potential V satisfies 



rju) G 



3M 



(4) 




1 Cjjmin{drixy^ , dr{x)~^} dr{x) 




if fi^O 
I'^l if M = 



(5) 
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for all X £ with dr{x) 7^ and all multi-indices 7] with jr/j < 3, and 

\V{x)-Zk\x-rk\-'\<Cr„^l + C (6) 

for |x— r^l < rmin/2 where rmin = niin^^^ l^fc— r^l- Note, in particular, that the Thomas-Fermi 
potential F^^(z,r, •) discussed in (35) below satisfies these requirements, by Theorem 20. 
So does the potential V{x) = — 1 (with M = 1, ro = 0, and dr{x) = \x\). 

The main new result in this section is the relativistic Scott correction to the semi-classical 
expansion for potentials of this form. It will be proved in Section 4 below. The power —3 in 
(5) is not optimal. 

Theorem 4 (Scott-corrected relativistic semi-classics). There exists a continuous, 
non-increasing function 5 : [0, 2/7r] — > M with S{0) = 1/4, such that for all h > 0, < P < 
, Tjs as in (3), and all potentials V : M.^ ^ M. satisfying (5) and (6) with rmin > ro > 
and maxjzi, . . . , zm} < 2/7r, we have 

„ M 
Tr [r^(-iW) - V{x)] _ - {2^hY^ / [i/ - V{v)\ _ dvdp - h'^ J2 zl^iP^^^h'^Zk) 

k=i 

< C/1-2+V10 . (7) 

Here, [x]- = min{x,0}. The constant C > depends only on M, ro, fj, and the other 
constants in (5) and (6). 

Moreover, we can find a density matrix 7, whose density satisfies (with \\ • Hg/s the 
L^/^ -norm) 



J p^{x)dx-2^/\3TrY^h-^ J \V{x)-f^dx 



< Ch-^+^/^ (8) 



and 

\\p^ - 2V2(37r2)-l/,-3|y_|3/2||^^^ < C^-2-1/10 ^ (9) 

such that 

„ M 

TT[iTf,{-ihV)-Vix)h] <(27r/i)-3 / [^^ - V{v)] _ dvdp + h-^^. ^l^(f^'^^^~' ^k) 

•' k=l 

+ C/i-2+Vio. (10) 

Remark 5. The term proportional to h~^ is called the Scott correction. If (3 = h^ then it 

only depends on the charges Zk, k = 1, . . . , M, of the Coulomb-singularities. Notice that the 
function in the semi-classical integral is the non-relativistic energy. This is also the reason 
why the leading Thomas-Fermi energy is independent of (3. 

Applying this theorem to the potential V{x) = — 1 (which satisfies (5) and (6) with 
M = 1, ro = 0, and dr{x) = \x\), and using a simple scaling argument, gives the following 
explicit characterization of the function S in Theorem 4 (see details in Lemma 27 in Section 4 
below) . 

Corollary 6 (Characterization of the Scott-correction S). The function S satisfies, 
uniformly for a G [0, 2/7r] , 

S{a) = \mJ^Tr:[Hc + k]_- {2tt)-^ j [Ip^ - \v\-'^ + k] _ dpdv^ , (11) 
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where 

/V-a-2A + a-4-a-2-|x|-i , a G (0, 2/7r] 

Remark 7. Another characterization of the function 5 is given in Lemma 25 in Section 4 
below. 

Remark 8. The result in Corollary 6 was proved in [24, Theorem 7.4], but only pointwise 
and only for a G (0, 2/7r). 



2. Preliminaries 

2.1. Analytic tools. We recall here the main analytic tools which we use throughout this 
paper. We do not prove all of them here but give the standard references. Various constants 
are denoted by the same letter C although its value may change from one line to the next. 
Let p > I, then a complex- valued function / (and only those will be considered here) is 

said to be in i7(R") if the norm ||/||p = (/ \f{x)\Pdxf^^ is finite. We denote by ( , ) the 
inner product on L'^iW"'); it is linear in the second and anti-linear in the first entry. For 
any 1 < p < t < q < oo wc have the inclusion f] L'^ C -L*, since by Holder's inequality 
ll/llt < ll/llpll/llg"'^ with Xp-^ + (1 - A)g-^ = t'^. We denote by / the Fourier transform 
of / G L2(M''), given by /(p) = (27r)-"/2 j g-ix py^^^ f^j. Schwartz functions on R", and 
extended by continuity to L'^{W^). 

We denote x_ = min{x, 0}, and let XA be the characteristic function of the set A; we 
write X = X(-oo,o] the characteristic function of (— oo,0]. We call 7 a density matrix on 
L'^{W^) if it is a trace class operator on L'^(W^) satisfying the operator inequality < 7 < 1. 
The density of a density matrix 7 is the L^-function such that Tr(70) = J p^{x)6{x)dx 
for all 9 G C^(M") considered as a multiplication operator. 

We also need an extension to many-particle states. Let tf) G 0^L2(R^ x {—1, 1}) be an 
iV-body wave-function. Its one-particle density is defined by 

N 

P^{x)='Y ^ ••• ^ \lp{xi,Si,...,XN,SN)\'^5{Xj -x)dXl■■■dXN■ 

j=lSl=±l sjv=±l 

The following two inequalities we recall are crucial in many of our estimates. They serve 
as replacements for the Lieb-Thirring inequality [22] used in the non-relativistic case. 

Theorem 9 (Daubechies inequality). One-body case: Letm > 0, f{u) = + Jti^ — 
m, and F{s) = /o*[/~^(i)]" c?i, where denotes the inverse function of f. Assume that 
V G Lj^^(R"), and let - A be the Laplacian in R". Then 

Tr[^/-A + m^-m + Vix)] _ > - C J F{\V{x)-\) dx , (13) 

where a;_ = min{x,0}, and C is some positive constant. 

Many-body case: Let V' £ A -^^^(R^ x {— Ij +1}) ^.^d let p^ = p be its one-particle density. 
Then 

N » 

(v., ^ [^-A, +m2 -m]i^)> G[p{x)] dx , (14) 
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where (with Cq = 0.326) 

G{p) = (3/8)m4Co<7[(p/Co)^/=^m-^] - mp , (15) 
with g{t) = t{\ + t2)i/2(i + 2*2) - log[t + (1 + ^2)1/2]. 

The asymptotic behaviour of G for smah, respectively large p is given by 

G{p) {2,/l0m)C,''''p^/\ G{p) ~ {3/A)C-"''p'lK (16) 

By a simple scaling, and using the definition of T and (16), respectively, we see that 

Tr [V-a-2 A + m2a-4 - ma^^ + v{x)\ _ (17) 

>-c»»/^/|vW_|.-/.,._c,,»/|v.(.,_|.«^. 

and 

N „ 

^ [y/-a-2Aj + m2a-4 - ma-^j > C mm{m-^ p{xf/^ , a-^pix)"^/^} dx . (18) 

i=i 

Both (17) and (18) also holds for a = 0, where wc let V— ci~2 A + m'^a~^ —ma~'^ = — A/2m, 
when a = 0. The original proofs of the inequalities (13) and (14) can be found in [6] (for 
a = 0, in [22]). 

Theorem 10 (Lieb-Yau inequality). Let n = 3. Let C > and R> and let 

2 



Hc,R = V^-^-C/R. (19) 



vr X 



Then, for any density matrix 7 and any Junction 6 with support in Br = {x\\x\ < R} we 
have that 

Tr[ejeHaR] > -4.4827 C^ii:'^{3/(47ri?3) j \e{x)\'^ dx} . (20) 

Note that when = 1 on Br then the term inside the brackets {} equals 1. 

We will need the following new operator inequality. The proof can be found in Appendix A. 

Theorem 11 (Critical Hydrogen inequality). Let n = 3. For any s G [0,1/2) there 
exists constants As,Bs > such that 

\^-^>^(-A)^-i?3. (21) 

We also use the following standard notation for the Coulomb energy, 

D{f) = D{f, f) = \jj(x)\x- y|- V(2/) dxdy . 

Theorem 12 (Hardy-Littlewood-Sobolev inequality). There exists a constant C such 
that 

D{f)<C\\f\\l„. (22) 

The sharp constant C has been found by Lieb [18]; see also [19]. It can be shown by 
Fourier transformation that / ^ ■\/ D(f) is a norm. This fact will play a role in the proof 
of the upper bound in our main Theorem 1. 

In order to localise the rclativistic kinetic energy we shall use the equivalent of the IMS- 
formula for the operator — A/2m. In the sequel, as before, \/— Q!~2 A + m'^a~^ — ma~^ = 
—A/2m, when a = 0. 
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Theorem 13 (Relativistic IMS formula). Let {9u)ueM be a family of positive bounded 
-functions on with bounded derivatives, and let dfi be a positive measure on M. such 
that Jj^ Ouixf dii{u) = 1 for all x G M^. Then for any f e H^/'^{R^), 

if, [V-a-^A + m^a-"^ - ma-^)f) (23) 

= [ {Ouf, (\/-a-2A + m2a-4 - ma-y^f) dfi{u) - (/, Lf) , 
Jm 

where the operator L is of the form 

[ Le^dti{u), (24) 
JM 

with Lg^ the bounded operator with kernel 

LeSx,v) = (27r)-2m2a-3|a; - ^1-2^2(^0-^ " y\)[Ou{x) - e.,{y)\\ (25) 

Here, K2 is a modified Bessel function of the second kind. For a = 0, Lq^ is multiplication 
by {Sldy)"^ jlm, where V— a^^A + rxv^ar'^ — moT''^ = —A/2m, when a = 0. 

A proof (and the definition of K2, and some of its properties) can be found in Appendix A. 
The following bound on the localisation error will be crucial. 

Theorem 14 (Localisation error). Let Q cM.^ and £ > 0. Let 6 be a Lipschitz continuous 

function satisfying < 9 < 1, dist(J7'^, supp V^) > i, and 9 is constant on $7*^. 

Then for all m > 0, a > there exists a positive operator Qg such that the following 
operator inequality holds: 

Le<Cm-^\\V9\\lxn + Qe, (26) 

with 

Tr[Qg] < Cma-^r^e-""'''^\\V9\\l^\n\ , (27) 

for a constant C > 0, independent of m,a,£,9, and O. Here, Xd and are the character- 
istic function and the volume, respectively, of the set J7. For a = 0, Qg = d. 

A proof can be found in Appendix A. Note that the first term, C m~^||V^||^X0) on the 
right side of (26) is similar to the error in the non-relativistic IMS formula for the operator 

— A/2m, except in this case one has ||V0||^XsuppV6i/2?Ti as the only error. 
When localising, we shall make use of the following. 

Theorem 15 (Partition of M"). Consider ip G Cq°(IR"') with support in the unit ball 
{\x\ < 1} and satisfying J ip{x)'^ dx = 1. Assume that £ : M"" — > M is a C^-map satisfying 
< £{u) < 1 and ||V^||oo < 1- Let J{x,u) be the Jacobian of the map u i— j^, i.e.. 



det 



{xj - Ui)dj£{u 
£{u) " 



J(x, u) = £{u) 

We set (fuix) = (p{j^) V J{x, u) ^('u)"/^ . Then, for all x G 



J 



ipu{xf £{u)-'' du = 1 , (28) 



and for all multi-indices rj G N" we have 
where Cr, depends only on rj. 



PV^Iloo < ^(u)-l''IC^ max PVIloo , (29) 
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This is Theorem 22 in [36]. 

We will consider potentials F : — M with Coulomb singularities of the form Zk\x — 
Rk\~^, k = 1, . . . , M, at points . . . , Rm € and with charges < zi,. . . , zm < 2/7r. 
Recall that (see (4); replace r by R) 

dn{x) = ram{\x-Rk\\k = l,...,M} , R = {Ri, . . . , Rm) e R^^ . (30) 

To treat such potentials we will need the following combination of Theorems 9 and 10. The 
proof can be found in Appendix A. 

Theorem 16 (Combined Daubechies-Lieb-Yau inequality). Let Ri, . . . ,Rm € R^, 

and assume W G Ljjjp(M^) satisfies 

W{x) > — - — 1-^ — Cvma'^ when d-R,{x) < am.~^, (31) 
dji{x) 

with au <2/ir and m > 0, a > 0, and dji as in (30). Assume also that the minimal distance 
between nuclei satisfies min^^^ \Rk — Ri\ > 2am~^. Then 

Tr [\/-a-2A + m2a-4 - ma"^ + W{x)] _ 

> -Cu^/'^a^/^m-Cm^/'^ J \W{x)-\^/^ dx - Ca^ J \W{x)-\Ux, (32) 

dji{x)>ani~^ d-[i{x)>am~^ 

where as before V— a^^A + m^a"^ — ma"^ = — A/2m, when a = 0. 

Finally, we come to the two inequalities which bound the many-body ground state energy 
in terms of a corresponding one-body energy. 

Theorem 17 (Correlation inequality). Let p -.M.^ be non-negative with D{p) < oo 
and let ^ : M.^ ^ R be a spherically symmetric, non-negative function with support in the 
unit ball such that f ^{x)dx < oo. For s > 0, let ^s{x) = s~^^{x/s). Then, for some 
constant C independent of N and s, we have^ 

N 

\xi-Xj\-^ >Y{p*\x\-^*^s)ixj)-D{p)-CNs-K (33) 
l<i<j<N j=l 

The proof can be found in Appendix A. 

Theorem 18 (Lieb's Variational Principle). Let ^ be a density matrix on L^(M^) sat- 
isfying 2Tr7 = 2 J p^{x)dx < Z (i.e., less than or equal to the total number of electrons) 
with kernel Pj{x) = j{x,x). Then 



E{Z,'R;a) < 2TT[{V^cP^^^A + or^ - a'^ - V{Z,K,x))-f] + D{2p^) . (34) 
The factors 2 above are due to the spin degeneracy, see [17]. 



^We denote convolution by *, i.e., {f*g){x) = f f{y)g{x — y) dy. We also abuse notation and write p* \x\ ^ 
instead of (p * | • \~^){x). 
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2.2. Thomas-Fermi theory. Consider z = {zi, . . . , zm) £ 1^+^ and r = {ri, . . . ,rM) S 
R^^. Let < p G L5/3(R3) n L^{R^) then the (non-relativistic) Thomas-Fermi (TF) energy 
functional, S'^^, is defined as 

£^^{P) = U^T^''?'^ j Pixf/'dx- J Viz,r,x)pix)dx + D{p), 

where V is as in (1). 

By the Hardy-Littlewood-Sobolev inequahty the Coulomb energy, D{p), is finite for func- 
tions p e L5/3(M3)nL^(M3) c L^/5(M3). Therefore, the TF-energy functional is well-defined. 
Here we only state without proof the properties about TF-theory which we use throughout 
the paper. The original proofs can be found in [21] and [16]. 

Theorem 19 (Thomas-Fermi minimizer). For all z = {zi,...,zm) £ 1^+^ and 
r = (ri, . . . , ta//) G there exists a unique non-negative p'^^{z,r,x) such that 

J p^^{z, r, x) dx = Yl!k=i o-nd 

fTF(^TF) ^ inf {fTF(^) I < p G L^/^{R^) n L\R^)} . 
We shall denote by E'^^{z,r) = £^^{p^^) the TF-energy. Moreover, let 

V'^^iz, r, x) = V{z, r, x) - p^^{z, r, •) * (35) 

be the TF-potential, then V'^^ > and p^^ > 0, and p^^ is the unique solution in L^/^(M^)n 
L^(R^) to the TF-equation: 

V'^^iz,r,x) = i(37r2)2/3p^^(z,r,x)2/^ . (36) 

Very crucial for a semi-classical approach is the scaling behavior of the TF-potential. It 
says that for any positive parameter h, 

V^^{z.r,x) = h^V^^{h-^z,hr,hx), (37) 

pT^(z,r,x) = h^p^'^{h-^z,hr,hx), (38) 

E^^{z,r) = h''E^^{h~\,hr). (39) 

By hr we mean that each coordinate is scaled by h, and likewise for h~'^z and hx. By 
the TF-equation (36), the equations (37) and (38) are obviously equivalent. Notice that 
the Coulomb-potential (the potential V in (1)) has the claimed scaling behavior. The rest 
follows from the uniqueness of the solution of the TF-energy functional. 

We shall now establish the crucial estimates that we need about the TF-potential. For 
each A; = 1, . . . , M we define the function 

Wfc(z,r,a;) = V'^^{z,r,x) - Zk\x - rk\~^ . (40) 

The function Wk can be continuously extended to x = r^. 

The first estimate in the next theorem is very similar to a corresponding estimate in [15] 
(recall that the function dr was defined in (4)). 

Theorem 20 (Estimate on V^^). Let z = {zi, . . . ,zm) G R+ and r = (ri,...,rM) € 
]^3M_ multi-indices r/ € N'^ and all x with dr{x) we have 

|ajyT^(z,r,x)| < C^min{dr(x)-\dr(x)-^}4(a;)-l''l , (41) 

where Crj > is a constant which depends on rj, zi, . . . , zm, and M . 
Moreover, for \x — rk\ < rniin/2, where rmin = rnin^^^ \rk — re\, we have 

-C<Wk{z,v,x)<Cr^^ + C, (42) 
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where the constants C > here depend on zi, . . . , zm, and M . 

Corollary 21 (Estimate on p^^ * * (^o — -^e^ $ : ^ M 6e a spherically 

symmetric, positive function with support in the unit hall and integral 1, and for t > Q, let 
$j(x) =i-3$(x/t). If p^^{x)= p^^{x,Y,x) then 

n < /,TF^|^|-i «TF < f CtininK(x)-V2^dr(x)-2} for dr{x)>2t 

0<p *\X\ -p *\X\ *$i<| ^^1/2 fordrix)<2t ^^^> 

with the function dr from (4), and some constant C > depending on zi, . . . , zm, and M. 

For the proof of (41) and (42) wc refer to [36]. (Note that in [36] it is claimed that 
Wk{z, r, x) > 0. This is not correct, but the proof in [36] does give that Wk{z, r,x) > — C.) 
The proof of (43) can be found in Appendix A. 

Remark 22. As is seen from the proofs in [36] and in Appendix A, the constants in 
Theorem 20 and Corollary 21 only depend on > when zi, . . . , € (0, zq\. 

The relation of Thomas-Fermi theory to semi-classical analysis is that the semi-classical 

density of a gas of non-interacting (non-rclativistic) electrons moving in the Thomas-Fermi 
potential 1/^^ is simply the Thomas- Fermi density. More precisely, the semi-classical approx- 
imation to the density of the projection onto the eigenspace corresponding to the negative 
eigenvalues of the Hamiltonian — — is 

\f Wr = ^'"i^-'r'iy^y^-^-^-) = p^"(->r,x). (44) 

ip2_v/TF(z,r,x)<0 

Here the factor two on the very left is due to the spin degeneracy. Similarly, the semi- 
classical approximation to the energy of the gas, i.e., to the sum of the negative eigenvalues 

of-iA-yTF^ is 

2 j [i/-F^^(z,r,.)]_|| = -l&j V^^iz,r,xf/Ux 

= E^^{z,r)+D{p^^{z,r,-)). (45) 

Since (by Theorem 20) the Thomas- Fermi potential V'^^{z,r, •) in (35) satisfies (5) and (6) 
(uniformly for zi, zm G (0, 2/7r]; see Remark 22), Theorem 4 implies that the density 
given in (44) and the energy given in (45) are the leading order terms also for the relativistic 
gas, i.e., for the operator T^(— i/iV) — V"^^ , < /3 < /i^, with as in (3). That the Thomas- 
Fermi energy is correct to leading order for r^2(— i/iV) — V"^^ was proved in [25]. Theorem 4 
establishes the first correction — the Scott correction — to the leading order. 

3. Proof of the relativistic Scott correction for the molecular ground 

STATE energy 

In this section we prove Theorem 1. Except for the correlation inequality we proceed 
in exactly the same manner as in the non-relativistic case [36]. In [36] correlations were 
controlled by the Lieb-Oxford inequality [20]. Applying this inequality, correlations can be 
estimated by the integral J p^l"^ involving the electronic density p. Using the non-relativistic 
Lieb-Thirring inequality such an integral can be seen to be of lower order than the total 
energy. In the present relativistic case the Daubechies inequality (14) a priori only allows 
us to conclude that the integral J p^l^ is of the same order as the total energy. We therefore 
follow a different strategy. 
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Proof of Theorem 1 (Lower bound). Let ip he a (normalised) ground state wave function 
and let s > 0. We will use the correlation inequality (33) with p{x) = p^^{Z,'R,x). Let 
be a function as in Theorem 17. Wc shall choose s = Z~^l'° . 

As above (see (4)) we have d-r{x) = min||a; — rk\ \ k = 1,...,M}. Note that for the 
physical positions of the nuclei we then have 

dR(x) =min{|x-i?fc| \k = l,...,M] = Z-^/^d^{Z^/'^x) . 

From the estimate in (43) with t = Z^/^s we obtain from the Thomas-Fermi scaling (38) 
that 

|p'^^(Z, R, ■) * \x\-^ - pTF(z, R, ■) * \x\-^ * $,(x)| < CZ^/'^s{g{x) + Z^l^) , 

where 

' (2.s)-V2 if di^{x) < 2s 
9{x) = I dR(a;)-V2 if 2s < d^{x) < Z-^l^ . 
^ if < d^{x) 

We find from the correlation estimate (Theorem 17) that 

(V^,//(Z,R;a)V'> 
z 

> J2 + - - R, xj)] ^) 

z 

+ <V', E iP^^C^' ^' •) * 1^1"' * ^«)(%)^> - D{p^^{Z, R, •)) - Cs-'Z 

i=i 

> 2Tr[V-a-2A + a-4 - a'^ - V^^{Z, R, x) - CZ^/^sg{x)] _ - D{p^^{Z, R, •)) 

- CsZ^/^ - Cs-^Z . (46) 

To control the error term with g above we shall use the combined Daubechies-Lieb-Yau 
inequality (Theorem 16) to estimate 



£Tr[V-a-2A + a-4 - a"^ - F^^(Z, R, x) - Ce'^ Z^l^sg{x)\ _ 

for some < e < 1 which we will choose to be e = Z^^/^. Wc use Theorem 16 with m = 1 
and f = max^. Z^. Then by assumption va < 2/7r. We must also check that the assumption 
(31) is satisfied, i.e., that for d^{x) < a wc have 

- V^^iZ, R, x) - Ce-^Z^''^sgix) > - — ^ - Cva'^ . 

dn{x) 

This follows from the definition of g and the estimate on the TF potential in (42) together 
with the Thomas-Fermi scaling (37) if 

^a<s< C{i^Z-^f{Za)-^e'^Z, r"?-^ + 1 < C{iyZ-^){Za)-^Z^/^, 

which, for Z large enough, is a consequence of the assumptions in the theorem and the 
choices of £ and s. Note, in particular, that vZ~^ = max/j z/j > (since J2k-^k = 1) 

and by assumption Za < minfc2/(7rzfc) < 2M/7r. The constants C above depend only on 
zi, . . . , zm, and M. 

According to the Thomas- Fermi estimate (41), the Thomas- Fermi scaling (37), the defi- 
nition of g, and the choices of s and e we have 

V^'^iZ, R, x) + Ce-^Z^/^sg{x) < Cmin{dR(a;)-^ Zdn{x)-^} . 



14 J. P. SOLOVEJ, T. 0STERGAARD S0RENSEN, AND W. L. SPITZER 

Thus the combined Daubechies-Lieb-Yau inequahty gives, since v < Z and Za < 2M/it, 
that 



£Tr[V-a-2A + a-4-a-2 - V'^^{Z,R,x) - Ce-^Z^/^sg{x)]- 

> -CeZ'^-Ce J {mm{dRixy^,Zdnix)-^}f^^dx 

-Cea^ j {m.ui{dn{x)''^,Zdn{x)~^})^dx 

> - Ce{Z^ + Z'^l'^) > - CeZ'^l^ . 

We return to the main term in (46). Using the Thomas- Fermi scaUng property (37) and 
replacing x by Z~^l^x we arrive at 



Tr [\/-a-2A + a-4 - q-^ _ Y^^ {7., R, x)\ _ 
= Z^/^^-^ TV[V-/?-^^2A + - - kV^^^\z, r,^)]_ , 
where we have chosen 

K = mm—> Za, h = i^'^Z-^/\ P = Z^'^a^n-\ (47) 

k TTZk 

We shall use /3 and h as the semi-classical parameters when we apply Theorem 4. It is 

therefore important that (3 < h^. This follows since p-^h^ = {Zay^K^ > 1. Note also that 
2/7r < n < 2M/7r since Zk < k = \, . . . , M, and Zk = 1- 

Putting this together with the estimate above into (46) we obtain (using the inequality 
Tr[X-|-y]_ > Tr[X]„ -|-Tr[y]_ for operators X and Y bounded from below (with a common 
core), and the choices of e and s) that 

> 2(1 - e)Z^/3^-i [^_/?-i/i2^ + ^-2 _ p-i _ ^v^Fi^z, r, x)] _ 
- CsZ'/^ - CsZ^/^ - Cs-'Z - D{p^^{Z, R, •)) 

> 2Z^/3k-i TV[V-/3-i/i2A + /3-2 - /j-i - KV'^^{z,r,x)]_ 
-L>(pTF(z,R,.))-CZ'-^/^ 

Now we apply the semi-classical result for potentials with Coulomb-like singularities from 

Theorem 4 to KV'^^{z,r, •) (recall that 2/7r < k < 2M/7r which ensures that the constants 
in (5) and (6) are uniform in k), and the calculation in (45). Then 



2Z^/\-^ Tr[V-/3-^/i2A + /?-2 _ p-^ _ kV^^{z, r, x)] _ 

M 

= ZV\E^\z,v) + D{p'^{z,v,.)))+2Y,ZlS{Zka) + 0{Z^-' 

k=l 

M 

= E^^{Z,R) + D{p^^{Z,R,.))+2j2zlSiZka)+0{Z'-'/''). 



k=l 



Note here that k cancels in the leading semi-classical term and in the Scott-term (the term 

with S). Also, 2/7r < K < 2M/tt ensures that the error is uniform in k. Here we have again 
used the TF scaling E^^^ {Z, R) = Z^/s^tf^^^ p) and D{p^^{Z, R, •)) = Zy^D{p^^{z, r, •)) • 
This finishes the proof of the lower bound. □ 
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Proof of Theorem 1 (Upper bound). The starting point now is Lieb's variational princi- 
ple, Theorem 18. By a simple rescaling the variational principle states that for any density 
matrix 7 on L^(M^) with 2 Tr7 < .Z we have 

£;(Z,R;a) < 2Z^/3 Tr[(V-a-2^-2A + a-4Z-8/3 - a-'^Z-'^'^ - F(z, r,f ))7] 

As for the lower bound we bring the TF-potential into play: 

Z-^/^E{Z,-R;a) < 2Tr[{V -a-^Z'^ A + a'^Z'^/^ - a-'^Z'^l^ - ^^^(z, r, f ))7] 

+ ZD{2Z-^p^ - p^^iz, r, •)) - ZD{p''^{-L, r, ■)) 
= 2k-i Tr [( V-^-i/i2A + /3-2 - - kV"^^{x, r, x))7] 

+ ZL>(2Z-V7-/''(z,r,-)) -^i?(/9^^(z,r,.)) , (48) 

where k, /i, and /? are chosen as in (47) in the proof of the lower bound. Note that with this 
choice of h and k we have from (36) that 

2^l\-i^h')-\KV^^{-L,v,x)f = V^(z,r,x)/2. 

We now choose a density matrix 7 according to Theorem 4 with V{x) = ^^'^^(z, r, x). 

Since / p'^^(z,r,x) dx = Efl^ = 1 we see from (8) that 2Tr7 < Z(l + CZ-^/3-i/i5) 
(recall that k"^ < 7r/2). Thus if we define 7 = (1 + CZ-V3-Vi5)-i^ see that the 
condition 2 Tr7 < Z is satisfied. 

Using the Hardy-Littlewood-Sobolev and (9) inequalities we find that 

ZD{2Z-'p~, - pTF(z, r, .)) < CZ-^\\p~, - Zp^^-., r, .)/2\\l^^ < CZV3-4/i5 ^ 
and thus 

Z£»(2Z-V7-/''(z,r,.)) < C(l + CZ-V3-i/i5)-2^^(2^-ip^_/F(^^^^.)) 

+ CZV3-2/15^(^TF(^^ ^ < C-Z^/^-VlS ^ (49) 

where we have used the triangle inequality for ^fD, and that D(p"'"^(z,r, •)) < C. 
Finally, if we use (10) and (45) we get as for the lower bound that 

2Z^/^K-^ Tr [(V-/3-i/i2A + /?-2 - - kT/TF (^^ r , x)) 7] 

M 

< E^^{Z,K) + D{p^^{Z,K,.))+2j2z!SiZ,a) + 0{Z'-'/''). 

i=i 

Since £^'^^(Z,R) > - CZ'^/^ and D{p^'^{Z,'R, •)) > we sec that the same estimate holds 
for 7 replaced by 7. This together with D{p^^ {Z,^,, ■)) = Zy'^D{p^^\z,r,-)) , (48), and 
(49) finishes the proof of the upper bound. □ 

The function S is continuous and non-increasing, and <S(0) = 1/4, according to Theorem 4. 
This finishes the proof of Theorem 1. 

4. RELATIVISTIC SEMI-CLASSICS FOR POTENTIALS WITH COULOMB-LIKE SINGULARITIES 

In this section we prove Theorem 4. The theorem will follow from using Theorem 23 below 
(a rescalcd version of the local semi-classical results for regular potentials in Theorem 32 in 
Section 5 below). We localise (Theorem 13) the operator using multi-scale analysis (The- 
orem 15), and control the localisation errors (Theorem 16). Near the singularities of the 
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potential, we compare with the Coulomb potential. To be able to do this, we first prove a 
Scott-corrected semi-classical result for a localised relativistic Hydrogen operator (Lemma 25 
below). The ingredients of the proof of the latter are the same (rescaled semi-classics, local- 
isation and multi-scale analysis, and estimating localisation errors). 

Theorem 23 (Rescaled semi-classics). Let n > 3 and let (p G Cq ^^(M") be supported in a 
ball Bi of radius £ > 0. Let V G C^{B() be a real potential, and let Tp{p) = I3~^p^ + /?~^ — 
I3~^ be the kinetic energy. Let Hj^ = T^(— i/iV) +V{x), h > 0, and (J(){v,q) = Tj^^q) + V{v). 
Then for all h,l3,f>0 with < 1, we have 



TV[0F^0]_ - (27r/i)^"y" (l){vfa(,{v,q)^dvdq 



^ n+6/5^n+4/5^n— 6/5 



where the constant C is independent of [5 and depends only on 

sup ||^l^l0''(/.||oo and sup \\f-^^\•^\^^V\\^ ■ 

H<n+i |r)|<3 

Moreover, there exists a density matrix 7 such that 

Tr[7(/)i/^0] < (27r/i)-" j 4>{vfap{v, q). dvdq + Ch-'^+^l^ f'^^^^t'"^'^ . 
The density p-y satisfies 

p^{x) - (27r/t)-"a;„|y_r/2(2 + p\V-\Y''^{x) < Ch-''+^/^^ f^'^/^^r^'^^ 
for (almost) all x G B(^, and 



(50) 
(51) 

(52) 
(53) 



(l){xfp^{x)dx-{2T^hyuJn I (^(a;)2|1/_|"/2(2 + /?|y_|)"/2(a;)dx 

^ ^j^— n+6/ 5 j;n— 6/5 £n— 6/5 



(54) 



where tOn is the volume of the unit ball Bi in M" . The constants C > in the above estimates 
again depend on the parameters as in (51). 

Proof. We introduce the unitary scaling operator {Utp){x) = Cr'"'/'^4){fr^x). Then 

U*4>[Tp{~ihV) + V{x)\4>U = f^i[Tf^p{-ihf-H-'V) + Vf4x)]ct>i , 
where (^^(x) = (j){ix), and Vf^£{x) = f~'^V{£x). Thus, 

Tr[4>Hp(l>]^ = f^Tr[(l>e[Tf^p{-ihf-H-^V) + Vf^x)]^ 

Note that (f)^ and Vf^e are supported in a ball of radius 1 and that for all multi-indices 77, 

||5V,||oo = ll^l^'aVlloo and ||W/,,|U=r'lKl''la^'V^||oo. 
Let af^e^p{u,q) = Tpp{q) + Vf^i{u). By Theorem 32 in Section 5 below there is a constant 
C depending on the parameters as in (51) so that, as long as < 1, 

Ti[cpH^<P]. - {27rhf-'r')-''f J Mufwiu,q)-dudq'^ < Cf^hf-'r')-''^''^' ■ (55) 
A simple change of variables gives 

(27r/t/-iri)-"/2 y Mufaf,i^piu,q)-dudq={27rh)-'' J ct>{vf a p{v , q) - dvdq , 
and we have proved (50). 
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Now, let 7//,/3 be the density matrix for 0^ [rgj2 (— i/i/^^^^^V) + Vf^i{x)](p£, which ac- 
cording to Lemma 34 satisfies 

f TV [0, [Tsp {-ihf-H-^V) + Vf,e{x)] (t^niAP 

< {2nhf-'r')-^f J Uuf<^f,e,^{u, q). dudq + C{hr^r^)-''+^"> , 

M^fp'yf,,Ax) dx-i2TThf-H-')-''^nJ Mx?\VfAx)-r^\'^+pf\yfAx)-\r^^dx 

The density matrix 7 = U^f^i^jsU*, whose density is P'y{x) = l~'^p^^ ^ ^{x/tj, then satisfies 
the properties in (52)-(54). □ 

Multi-scale Analysis. The rescaled semi-classics of Theorem 23 will be used in balls of 
varying radius. This idea goes back to Ivrii [15, 14]. We introduce a variable scale ^ and a 
corresponding family of localisation functions {^Pu\u&?i which will also be used in the proof 
of Theorem 4. 

Definition 24 (Scale for multi-scale analysis). Let < £0 < 1 be a parameter that we 
shall choose explicitly below, and let ri,. . . ,rM £ Define 

£(x) = i(l + ^(|x-r,|^+£grV2) . (56) 

k=l 

Note that ^ is a smooth function (due to ^0) with 

0<£(x)<l/2 and ||V£||oo < 1/2 . (57) 
Note also that in terms of the function d = dr from (4) we have 

1(1 + UyHo < 1(1 + M{dixf + ilr^^Y^ < iix) < ^{d{xf + ilf^ . (58) 
In particular, we have 

^(x) > i(l + M)-imin{(i(x),l}. (59) 

We fix a localisation function ip G Cq°(]R^) with support in {\x\ < 1} and such that 
J" ip{x)'^ dx = 1. According to Theorem 15 we can find a corresponding family of func- 
tions (fu G C^(M^), u € M^, where (p„ is supported in the ball {\x — u\ < £{u)}, with the 
properties that 



/ 



ipuixfi{u)-^ du = 1 and ||9V«l|oo < C^(«)"l''l , (60) 



for all multi-indices r], where C > depends only on ij and tp. For d{u) > 21q we have i{u) < 
^/5d{u)/4 and hence for all x with jx — u| < £{u) we have (note that d{u) < d{x) + \x — u\ 
and a/5/4 < 1) that 

e{u) < d{u) and d{u) < Cd{x) . (61) 

As a first step towards the Scott correction for Coulomb-type potentials we deal with 
the Scott correction for a single relativistic Hydrogen atom. This method for proving the 
existence of a Scott correction in the semi-classical expansion of the sum of eigenvalues 
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of an operator with a (homogeneous) singular potential without explicitly knowing the 
eigenvalues was first used by Sobolev [33] when studying (non-relativistic) operators with 
magnetic fields. 

Lemma 25 (Scott-corrected localised Hydrogen). There exists a non-increasing func- 
tion S : [0,2/7r] M, with 5(0) = 1/4, such that, if (j)r{x) = 0(a;/r), r € (0,oo), with 
G C^(M3), Q<<j)<l, satisfying y/T^ G C^iR) and 

1 for \x\ < 1 
for \x\>2 ' 



then there exists C > depending only on cj) such that for all a G [0, 2/7r] and r G (0, oo), 



TY[^rHc{a)^r]- - (27r)-3 j ^vf - \v\-^] _ dvdq - S{a) 



< Cr-^/^° , (62) 



where 



JV-a-2A + a-4-a-2-|x|-i , a G (0, 2/7r] 

As emphasised in Remark 5, the function in the semi-classical integral in (62) is the non- 
relativistic energy. See also Lemma 27 below for an alternative description of the function 
S. 

Remark 26. A result similar to the one in Lemma 25 was proved in [24, Theorem 7.1], but 
without uniform control in a and only for a G (0, 2/7r). 

Proof of Lemma 25. We fix a G [0, 2/7r] and write He = Hc{a). We define for r > 

Sr = TT[(^rHcia)(t>r]_ - i'^^)'^ J (t>rivf[h'^ " \v\~'^] _ dvdq . (64) 

We will show that Sr has a limit as r ^ oo. 

Let R > 2r. We estimate the difference between Ti[(PjiHc4>r]- and TT[(prHc4>r]- semi- 
classically. The leading semi-classical term involves the relativistic energy which is then 
compared to the non-relativistic energy. Below all constants will depend only on and in 
particular not on a G [0, 2/ it]. 

Denote V'r = y^l — (pr- t^*^ relativistic IMS formula (23), 

He = (prHcCpr + i^rHclpr - - L^^ , 

where L^^ and L^^ are given by (24) and (25) {M = {1, 2}). We multiply with (pR and get 
that 

(f)RHc(f)R = 4>rHc(t>r + (t>RAHcA(f>R - 4>R{L'<t>r + L^r)4>R ■ 

We have used that (t>R4>r = (pr since R > 2r. Now, let jr = x{^rHc^r) be the projection 
onto the negative part of (^rHqc^r. Then, by the variational principle and Theorem 14 (with 
m = 1, £ = r, = i?(0, 3r), and 6 = (pr and ipr, respectively), 

Tr[(PRHc(pR]- 

= TT[yR(j)rHc(l)r] + Tr[yR(pRiprHcipr4>R] - Tic[yR(pR(L^^ + L^r)4>R] 
> TT[-fRMHc-Cr-^)A] + T^y^bR^RMHc-Cr-^hr)i^rM (65) 
- Cr-^ . 
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Here, C is independent of a. We treat the part of the locaUsation error coming from the 
first term in (65). We spHt Hq = (1 — e)Hc + eHq for some < £ < 1 to be chosen and use 
the second term to control the error term. 

By Theorem 16 (with M = 1, i?o = 0, d{x) = \x\, m = 1 and v = 1), 

a<\x\<2r a<|s|<2r- 

assuming e^^r^^ < Ca^^ and using that a < 2/7r. We may choose e = r^^ if we assume 
that r > 1 (note that then indeed s~^r~^ = r"^ < 1 < 2a~^/7r). We then obtain 

As a result, we have shown that 
Tr[(j)RHcM- 

> (1 - e) Tr[7^ 0.Fc0r] + Ti[-fR^RMHc - Cr-^MA^n] - Cr-^'^ 

> TT[cbrHc(t>r]- + Tv[^rA{Hc " Cr-^ ct>^r)AM- " Cr~^^^ . 

We will treat the term Tr[(pR'tpr{Hc — Cr~'^(f)3r)ipr(pR]- by our semi-classical estimates 
in Section 5 below. We first rescale. Define the unitary scaling operator {Uip){x) = 
R-^/^(p{R-^x). Then 

He : = U*{Hc-Cr-^<p3r)U 

= i?-^(\/-Q-2A + R^a-'^ - Ra-^ - \x\-^ - C Rr'^ 4)^^ / r{x)) 

= R-\Tp{-\hV) - \x\-^ - CRr~^4>^r/R{x)) (66) 

with (3 = a^R^^ (< R^^) and h = R'^/"^. Let (j)R^r = (t>R'4'r = 4>rV^ ~ '^r ^(^) = 
4>R,r{Rx) (see (3) for Tp). In this way, (pRipriHc — Cr~'^(l)sr)ipr4'R V'-^^cV' unitarily 
equivalent. 

Now, let i and be the functions in (56) and (60), respectively, when M = 1, ri = 0, 
and £o = h'^ = R~^. By another relativistic IMS localisation we get that 

i;Hci^ = R-^ J ij(pu{Ti3i-ihV)-\x\-^ -CRr-^(t>3r/Rix))vui^eiu)-^du 

r/3R<\u\<5/2 
r/3fl<|u|<5/2 

We have used that Vv'u = for \u\ [r/3i?,5/2], since l{u) < ^d^P + ^o)^^^ (see (58)) and 
suppV' C {r/R < \x\ < 2}, swppipu C {|a; — n| < £{u)}. 

Concerning L^^, Theorem 14 with £ = £{u)/2, m = R, and $7 = r2„ = — -ul < 3£(u)/2} 
gives 

L^^<CR-'£{u)-\n^ + Q^^, 
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with 

Here we have used (60) . 

Notice that if the supports of tpu and tpui overlap then [-u — n'| < £{u) + £{u') and thus 

i{u') < £{u) + ||V^||oo(^M + i{u')) . (68) 

Therefore, since ||V£||oo < 1/2, we have that e{u') < C£{u) and thus l{u)-^ < C^(n')"^ 
Similarly, i{u) < Ci{u'), and so Xfiu ^ X{\x-u\<ce{u')} if the supports of ipu and ipu' overlap. 
Using this and (60) we get for all x G 



\-3 



du 



J {i{u)-\nAx))£{ur'du = j {£{u)-\nSx)){j vl{x)l{u')-^du')^{u)- 
< C j Vu'{x)l{u')-'^^u'{x)£{v!)~^ dv! . (69) 
Rewriting the last integral with u as integration variable we get 

-R~^ j i^Q,p^ipi{u)-^du. 

Here we have also used that Rr~'^ (p^j. / ji{x) < Ch'^i{u)^'^ for x in the support of This 
is a consequence of £{u) < ^\u\ + ^£q < ^\x\ + ^£{u) + £o for x in the support of ipu which 
implies that £{u) < \x\ + io < Cr/R for x in the support of aiid (f)3r/R- 

We will now use Theorem 23 (with cf) = ipfu, £ = £{u), = {\x — u\< £{u)}, f = f{u) = 
\u\~^/'^) on 

ijiPu{Tp{-ihV) - \x\-^ - Ch''£iu)-^)ipui^ , 
for u fixed with \u\ G [r/3i?, 5/2]. We claim that 

\m{^iPu)\\oo < C/(u)-l''l for all r]eN\ (70) 

This follows from (60), (61), and the estimate \d^'4>{x)\ < C^|x|~l^l. It suffices to check 
the latter for 1 < |a;| < 2 and r/R < \x\ < 2r/R, due to the support properties of ij:. 
Furthermore, for r > 3, \x\^^ + Ch?t{u)~'^ is smooth (as a function of x) on (use (58), 
£o = R~^, and \u\ > r/3R), and satisfies 

sup |92(|xr^ + Ch'^£{u)-'^)\ < Cr^f{uf£{u)-\'^\ for all G N^ (71) 

\x-u\<t{u) 

with f{u) = \u\'~^^'^ . For the Coulomb potential, this is trivial. For the other term, only the 
statement for r/ = is non-trivial; it follows from (59), h = R^^^"^, and \u\ > r/3R. Finally, 
the condition j\uf'l3 < 1 is also satisfied (when r > 3), since \u\ > r/3R and /? < R~^. 
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From Theorem 23 we conclude that 



> R-\2Trh)-^ j i^{vf^u{vf [Tp{q) - \v\-^ - Ch^i{u)-^]- i(u)-^ dudvdq 

r/3R<\u\<5/2 

-CR-'h-^+'/' J f{u)^-'/'i{u)-'-'/'du 

r/3R<\u\<5/2 

- R-^ J Tr[tPQ^„tP]i{u)-^du. 



r/3R<\u\<5/2 

Integrating the semi-classical error using f{u) = \u\~^/^, (59), and R > r gives the lower 
bound -Ci?-i/i"2+i/5(^/^)i/io = _ Cr-i/^o. 

From (67) it follows, using (59), a < 2/tt, and R> r, that 



R-^ J TiliPQ^^iP] £{u)-^ du < C J a-^e-"" '^^''y^i{uy 



''du 



r/3R<\u\<5/2 



r/3R<\u\<5/2 

< Cr-^e-^/^ 



Since supptpu C {t; | |ti — ?;| < i{u)} and |m| < 5/2 we have \v\ < \u\ + £{u) < Ci{u) on 
supp (fu- Using this, integrating in u (using (60)), we get 



R-\2nh)-^ J jp{vfipu{vf [Tp{q) - \v\-^ - Ch^iiu)-^]- i{u)-^ dudvdq 

r/3R<\u\<5/2 

> R'^J^j^ j ^i^f WP'^Q^ + - - - Ch^\v\-^]- dvdq . 

In order to compare this latter relativistic semi-classical expression with the non- 
relativistic semi-classical one we use the inequality — y_| < |a;— y| and a Taylor expansion 
of + 1 — 1 to arrive at 

/ I - a] _ - Wp-^q'' + a - 6]_ I dq 

< CP{(3{a + bf + 2(a + h)f'^ + Ch{f3{a + hf + 2(a + h)fl'^ (72) 
for all a, 6 > 0. This gives, using h?' = R~^ and (3 < R~^, that 

J ^P{vf _ \v\-^]_- [y/p-^q'^ + /3-2 - - \v\-^ -Ch^\v\-^]_^ dvdq 

< CR-^ j \v\-''l'^dv < C{Rr)-^/^, (73) 



r/R<\v\<2 



since R> r > 1. 

Thus undoing the scaling we arrive at 



R-^{2Trh)-^ J ipivfipuivf [Tp{q) - \v\-^ - Ch^e{u)-^]^ e{u)-^ dudvdq 

2 

> (27r)-3 J cj>R,r{vf [\q^ - W^]- dvdq - Ct-^'^^ . 



r/3R<\u\<b/2 



22 J. P. SOLOVEJ, T. 0STERGAARD S0RENSEN, AND W. L. SPITZER 

Summarizing, we have proved that there exists a constant C = C {(/)), independent of a G 
[0, 2/7r], such that for r large enough, and R > 2r, 



> TiicPrHccPr]- + (27r)-3 J cPn^riv)^ _ dvdq - Cr-^'^"^ . (74) 

Next, we want to bound Tr[0ijiJc<?^ii]- from above by Tr[^r-ffc^r]- by constructing a 
density matrix. To this end, we first set 7,. = x{4'rHc4'r)- Then we let 7„ be the den- 
sity matrix we get when we use Theorem 23 for the rescalcd operator tpipuHQ(puip (now 
with He = U*HcU with U as in (66)), for fixed u with \u\ G [r/3i?, 5/2], and set 
7u = U(pulu'PuU*. Finally, we define 

7 = (/)r7r^r + J i'rlu i>r i{u)~^ du . (75) 
r-/3i?<|u|<5/2 

Since < 7 < 1 and / ipl{x)l{u)-^ du = 1, 



< y ju^iuy^du < 1, 



and so we see, by multiplication with ijjj. on both sides, that < 7 < 1. Also, 7 is trivially 
trace class. By the variational principle we obtain that 

Tr[(f)RHc(l)R]- < Tr[^RHc(f)m] 

+ j TT[i;r<t>RHc(t>R'4^rlu]^{u)-'^ du 
r/3R<|ii|<5/2 



r/3R<|u|<5/2 

Here we have used that (t)R(t)rHc(t>r<i>R = <i>rHc<i>ri since R > 2r, and again scaled the 
operators inside the trace in the integrand. Using Theorem 23 we can bound the integral 
from above by 

R~^{2Trh)-^ J ^pivf^uivf[T0iq) - \v\-^] J{u)-^ dudvdq 

+ ci?-i/i-2+v5 J f{u)^-y^iiu)-'-y'du. 

r/3R<\u\<5/2 

As in the case of the lower bound, the error term is bounded by Cr~^/^°. 

Integrating with respect to u in the semi-classical expression above, changing back coor- 
dinates, and using (73), we conclude that 

Tr[(l)RHcM- < Tr[cl>rHc(l)r]- + (27r)-^ j (j)R,r{v? [h^ - 1^1"'] _ dvdq + Cr-^'^'' . (76) 
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Combining (74) and (76) we have shown that for R > 2r, 
\Sr — Sr\ 

< \Tr[4>RHc(PR]- - 'IV[0,ifc</)r]- + (27r)-3 J {Mvf - M^f) [W - 1^1"'] - dvdq 

Hence, {iSnlneN is a Cauchy-sequence and with S = S{a) the hmiting value we have 

This proves (62). That S is non-increasing follows from the fact that T^2(p) (see (3)) is 
decreasing in a. Finally, that 5(0) = 1/4 is a well-known fact [36]. □ 

Proof of Theorem 4- Using the combined Daubechies-Lieb-Yau inequality (see Theorem 16) 
with a = /3^/^/i~^(< 1) and m = h~'^ we may assume that h is bounded by some constant, 
which we may choose small depending on M and ro, using that Zk < 2/7r, k = 1, . . . , M, 
and that 5 is a bounded function (since it is non-increasing; see Lemma 25). 

In order to control the region close to and far away from all the nuclei we introduce 
localisation functions 9± G C^(M) with the properties that < 9± < 1 and 

(1) e^_ + el = 1, 

(2) e^{t) = liit<l and ^_(t) = for t > 2. 

Let < r < ro/4 and < ro < -R and define ^±{x) = 9±{d{x)/R) and ^±(x) = 0±{d{x)/r) 
(with d = dr as in (4) ) . We choose (assuming h is small enough) 

. = rV and «={^f- , (77) 

where 6 = h < 1/2 and i?^ = C/x~^ is chosen such that — V{x) > for d{x) > (see 
(5)). We will keep writing 5 and R in the calculations below to show why these choices are 
optimal. Clearly, ^'i_ + ^\ = 1, -|- 0^ = 1, and (fPi + + ^\ = 1. Note also that 

M 

<t^-{^) = ^^r,k{x) with er,k{x) = 9-{\x-rk\/r). 
k=l 

Step 1: Lower bound on the quantum energy. 

By the relativistic IMS formula (23) and Theorem 14 with m = h^^, a = P^^'^h~^{< 1), 
and either i = R, Q = {d(x) < 3i?}, and 9 = $± respectively, or i = r, O, = {d(x) < 3r}, 
and 9 = 9r,k, k = 1,. . . , M, or 9 = respectively, we find that 

Tpi-ihV) - V{x) 

= ^+{Tf){-ihV) - V{x))^+ + $_(T^(-iW) - F(x))$_ - - 

M 

= Y.'^r,k{Tp{-ihV) - V{x))9r,k + $-0+(r^(-iW) - F(x))0+$_ 
fc=l 

M 



+ (r^(-iW) - F(x))$+ -^-{Y. + ^4>+)^- - - ' (^^8) 



fc=l 



with 



< C/i^||V$±||^XKx)<3K}+Q*±, (79) 
Tr[g$±] < CI3-^R-^e-^^"'^'^~''^\\V^±\\l,\{d{x)<m}\, (80) 
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and (with, by abuse of notation, = X^^i ^9rk) 

L<f>^ < Ch'^\\V(t>±\\l^X{d{x)<3r} +Q<t>±, (81) 
Tr[Q^J < Cp-'r-'e-^P'''''^~''\\V<P±\\l\{d{x)<3r}\. (82) 

Using |{(i(.T) < 3R}\ < SGttMR^, ||V^>±||oo < CR~^ (and the corresponding estimates for r 
and (f)±), P < h"^, and h small, it follows that 

Tr[g$±] < Ch^R-^e-^-'"^'^ < CnH^ , Tr[Q^J < Ch\-^e-''~''/^ < CnH'' , 

for any > by the choices (77). 
Hence we have that 

Tr[T^{-ihV)-Vix)]_ 

M 

> ^Tr[e,,fc(r^(-iW) - V{x) - ch\-^)er,k]_ 

k=l 

+ iv[$_,^+(r;3(-iw) - v{x) - ch\-\{d{x)<3r} - ch^R-^)(l^+^-]_ 

+ Tr[^+{Tf,{-ihV) - V{x) - Ch^R-\^ai.)<3R})^+] _ - C^-^+Vio . (33) 

Each of the first three terms above will be compared to the corresponding semi-classical 
expression. We shall treat the three terms by different methods. The first term will be 
calculated using the Scott correction for Hydrogen in Lemma 25. The second term will be 
computed using the local rescaled semi-classics in Theorem 23. The last term is an error 
term which wc will treat first. 
Control of the third term in (83). 

We use the Daubechies inequality (17) with m = and a = /3^/^/i~^(< 1). In the case 
H = we obtain, using the choice (77) of R, 

TV[$+(r^(-iW) - V{x) - Ch''R-\{ai.)<3R})^+]^ 

> -Ch-^M [ \x\-^^/^dx-CM I \x\-^^dx-Ch^R-^ -Ch^R-^ 



\x\>R \x\>R 

> - C{h-^R-^I^ + + h'^R-^ ~ h^R-'>) > - Ch^/^ . (84) 

The case n ^ gives a smaller error since — y > on the support of $_|_ in this case. 
Control of the first term in (83). 
Using (6) and (77) we have 

M 

^Tr[er,k{Tp{-ihV) - V{x) - Ch\-yr,k]- 
k=l 

M 

> Y,^[^r{Tpi-ihV) - Zk\x\-^ - CS'^h-^)er]_ , 

k=l 

where we have written 6r{x) = 9^{\x\/r). We have used here that 

Cr-l + C < Crg-i + C < C5^h-^ . (85) 

It is this relation which sets a lower bound on 6. We will control the error using the combined 
Daubechies-Lieb-Yau inequality in Theorem 16 with m = h^^ and a = j3^^'^h~^{< 1). Note 
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that ma~ -1 = /?-V2/j-i > h-2_ Thus using Theorem 16 we find, for ah density matrices 7 
and all s > S'^, that 

eTr [-f{er{Tf3{-ihV) - Zk\x\-^ - Ce-^5'^h-'^)9r)] > - CieS'^/^ + e'^^^S^ + e-^S^)h-^. 
Thus for all density matrices 7 and all e > 5^ we have 

Cd'^h-^Ti^ye^] < £Tr [jer{Tp{-ihV) - Zk\x\-^)er]+C{ed-^/^ +e-'^/^6'^ +e-'^6^)h-^ . (86) 
Hence 

M 

J2H^r,k{TpHhV) - Vix) - Ch\-yr,k]- (87) 

k=l 

M 

> {1 - e)^Tr[er{Tp{-ihV) - Zk\x\-^)er]_ - C{ed-^/^ + e-^/H^ + e-H^)h-\ 

k=l 

For the corresponding semi-classical expression we have from (6) and (85) (using S < 1/2 
and \x- — y-\ < \x — y\) that 

M 



{2Trh)-^J (t)-{vf [\p^ - V{v)\ _ dvdp - ^{2'Kh)-^j Brivf [\p^ - Zk\v\-^] _ dvdp 

k=l 

A simple rescaling applied to the local Hydrogen result in Lemma 25 gives that 
Tr[0^(T^(-iW) - Zk\x\-^)9r]_- {2TTh)-^j Orivfllp^ - Zk\v\-^] _dvdp 

-zlh-^S{p'/^h-'zk) < Ch-\h-\)-'/'^ = (89) 

Combining (87) , (88) , and (89) , using that 5 is a bounded function (since it is non- increasing; 
see Lemma 25), that S < 1/2, and that 

(27r/i)-3 J er{vf [\p^ - Zk\v\-^] _ dvdp < Ch-^r^/^ = Ch-^S-^/^ , 

and choosing e = 5, we conclude that 

Tr[0_(r^(-iW)-y(f)-C/iV-2)(^_]_ > {2'Kh)-^j (l)-{vf[\q^ -V{v)\_dvdq 

M 

+ h-^ 4S{(5^/^h-\k) - C<5Vio^-2 (90) 

k=l 

Control of the second term in (83). 

Here we use the local rescaled semi-classics in Theorem 23. Before we apply our semi- 
classical estimates on the support of $-0+ we localise using the functions tpu from (60) for 

general M and with l{u) as in (56), with = r/4. From (77) and the choice of 8 it follows 
that £o < 1 for h small enough. If x is in the support of $-^+ and in the support of (pu 
then d{u) > r/2 = 2^o since (using (58)) 

r < d(x) < d{u) + £{u) < d{u) + max{d(M), £o} , 
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and also d{u) < 2R + 1 since £{u) < 1/2. Using again the relativistic IMS localisation (23) 
we thus have 

(r;3(-iW) - v{x) - ch\-''x{d(x)<3r} - ch'R-'')(t>+^- 

= J $_0+<^4r^(-iW) - Vix) - C/iV-2x{d(,)<3.} - Ch^R-^)vu(t>+^- iiuy^ du 

r/2<d{u)<2R+l 

- j^_4>+L^^4>+^_l{ur^du. (91) 

r/2<d(u)<2R+l 

Concerning L^p^, Theorem 14 with ^ = i{u)/2 and = = {\x — u\ < 3£(u)/2} (and 
m = h~'^ and a = 0^/'^h~^{< 1)) gives, using (29), that 

L^^ < ChH{ur\a^ + , (92) 

with 

nQ^u\ < C/3-^e-(^^'''^)-^^H < Cat^^ (93) 
for ah A?^ > as a consequence of i{u) > (1 + M)~V/8 and /? < h"^. Thus 

< Cjv/i^ for all iV > . 



Tr 



r/2<d{u)<2R+l 

By the same arguments as in the proof of Lemma 25 above (see (68) and (69)) we join 
the new localisation error term (from (91), (92)) with the previous localisation errors from 
(79) and (81). Since i{u) < max{d{u) , r / 4} we have R''^ < Ci{u)-^ for d{u) <2R + l (and 
h small enough when /x = 0; for /x 7^ 0, use £{u) < 1/2) and, by (61) (valid on the support 
of ipu when d{u) > r/2 = 2£q), 

r~'^X{d(x)<3r}{x)fuixf < Ce{uy'^<f.^{xf . 

This way, we have proved that 

Tr (r^(-iW) - V{x) - C/iV-2x{d(,)<3r} - Ch^R-^)^+^-]_ (94) 
> y Tr [(P+^u {Tp{-ihV) - V{x) - Ch^i{u)-^) (pu4>+] _ iiu)'^ du - C/j-^+Vio _ 

r/2<d{u)<2R+l 

Note that there is no need to write on the right side, since in general Tr($^$)_ > Tr^_ 
for any self-adjoint operator A and any function < $ < 1. 

For u such that d{u) > r/2 = 2io and d{u) < 2R + 1 we have from (5) and (77) that 

sup |9''(y(x) - Ch'^e{u)-'^)\ < Cf{ufi{u)-\'^\ for |r?| < 3 , 

\x-u\<i{u) 

W{<t>+<t>u)\\oo<Cr,l{u)-\'^\ for |r/|<7, 

with 



^^""^ ~ [ min{d(u)-V2,d(u)-3/2} if ^ = • ^^^^ 

We have also used that d{u) > b'^h^ 12 > h? and min{l, (i(n)} < Cl{u). 

We are therefore in a position to use the rcscaled semi-classics in Theorem 23 on the ball 
{\x — u\ < f\ with £ = £(ti), / = /(ti), and (p = <p+<pu for each u with r/2 < d{u) < 2R + 1. 
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Note in particular that Pf{u) < (5d{uY^ < 2(3/r = 2/3()/i"^ < 25 < 1 . We conclude that 
for ah u with r/2 < d{u) <2R + 1, 

lh:[<t>+^u{Tp{-ihV) - V{x) - Ch't{ur^)^u<t^+\_ 

- {2'Kh)-^ J (f)+{vfipu{vf [ V/3-ig2 + /3-2 _ ^-1 _ v{v) - Ch'^i{u)-^] _ dvdq 
< C^-2+V5/(n)4-V5^(y)2-l/5 _ (9g) 

The semi-classical integral may be estimated using (72) 

J [Vf3-W + - r ' - V{v) - ChHiu)-^] ^dq- J [y^-Viv)]^dq 

< Ch^{\V{v)\ + /i2£(^;)-2)V2 ^ chH{v)-^{\V{v)\ + h^e{v)-^f^^ 

< C{h'^\V{v)\''/^ + h'^e{v)-'^\V{v)f/^ + {M{v)-^f) , (97) 

for V in the support of ipu, since then wc have £{v) < 3£{u)/2 (see (57)) and |l^(i')| < 
Cdivy^ < Cd{u)-^ < Ch-'^ (see (5), (61) and (58)). We have also used that (3 < h"^ and 
that, by (58) and (77) {Iq = r/4), h?l{v)-^ < CK^r-^ = CS'^h-'^ < Ch-^. 

Combining (94), (96), and (97) (remembering that d{u) < Cd{v) if v is in the support of 
ifu and d{u) > r/2 = 2£o) we find, using (60), (5), and (95), that 

TV[$_0+(r;j(-iW) - V{x) - Ch\-\{di.)<3r} - Ch''R-^)4>+<^.]_ 

> {2iTh)-^ j (l)^{vf[\q^ -V{v)\_dqdv-C j h'^+^Z^ f{uy^/H{u)~^/^ du 

r/2<d{u)<2R+l 

-C J h-^h'^d{v)-'^/^ + h^£{v)-'^f{vf + {M{v)-^f)dv-Ch-^+^Z^^ . (98) 

C-''r<d{v)<2R+2 

If /X 7^ the error term in (98) is controlled as follows: 

I {h-^+^"'f{vf'''H{v)-^"> + h-H{v)-"^ + h-H{v)-^f{vf + hH{v)-'') dv 

C-^r<d(v)<2R+2 

< C J (/l-2+V5|^,|-19/10^in|l,|^,|}-6/5+;,-l|^|-7/2 
C-'^r<\v\<2R+2 

+ min{l, |-y|}"^h^r^/^ + h'^ min{l, |?^|}"^ ) dv 

< C(/l-2+V5(i?lVlO + ^-1/10) + /,-V-l/2 + ^-1^3/2 ^ ^2^3 + ^2^-2) 

< (99) 

with the choices (77) where i? = i?^;^ is a constant. 
If /X = we get instead 



/ 



(/l-2+V5/(^;)19/5^(^)-6/5 ^ ^-l^(^)-7/2 ^ ^-l^^^^-2^(^)3 ^ ;,2^(^)-5) 
C-^r<d(v)<2R+2 

< C(/t-2+V5^-l/10 + /,-V-l/2 + ^2^3 + ^2^-2) < Ch-^+inO_ (^qq) 
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If we insert the last two estimates (100) and (99) into (98) and then together with (84) 
and (90) into (83) we arrive at a lower bound on the quantum energy corresponding to one 
direction in (7). 

Step 2: Upper bound on the quantum energy. 

We obtain an upper bound on the quantum energy by choosing the density matrix 



M 

7 



M 

^Or,klk&r,k + 1 <l>+^ulu^u<i>+ ^{u)~^ du , (101) 



d{u)<2R+l 



where 7^, k = 1, . . . , M, are the density matrices 

Ik = x{dr,k{Tf){-ihV) - Zk\x - rk\~^)er,k) 

and 7tj, u G M^, are the density matrices given in Theorem 23 for the potential V with B£ 
being the ball centered at u, i = £{u), f = f{u) (see (95)), and (p = (p+(pu- Since 

M 

^e2, + 0^ = </>2_+</.2 =1, (102) 
k=l 

we immediately see from (60) that 7 is a density matrix. 

Using this density matrix as a trial state we obtain from Theorem 23 that 

Tr[r^(-iW)-y(x)]_ 

M 

< J]iv[0,,fe(r^(-iW)-F(f))0,,fe]_ 

k=l 

+ (27r/i)-^ J (f)+{vf(puivf [V^-W + - - V{v)] Jiu)-^ dvdqdu 

d{u)<2R+l 

+ C^-2+V5 J f(uf^l^t(u)-^/^du, (103) 



r/2«i{u)<2R+l 



where we have used the fact that 0+ and have overlapping supports only if d{u) > r/2. 
The last error term is estimated by Ch~^'^^^^^ as in the lower bound. 

Using that + /3~^ — < and the normalization of (pu (60) we find that 

(27r/i)-3 j <l)+{vf(pu{vf [VP~^Q^ + - - yiv)] J{u)~^ dvdqdu 

d{u)<2R+l 

< j (l)^{vf[\q^ -V{v)\_dvdq 

- (27r/i)-3 j (t)+{vf^u{vf [\q^ - V{v)\ Jiu)"^ dvdqdu 



d(u)>2R+l 

< {2T:h)-^ j <t>+{yf[\q^ -V{v)\_dvdq + Ch-^ j \V{v)-\^l^ dv . 

d{v)>2R 
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If / the last term vanishes by the choice oi R = R/^. If /Li = it may be estimated using 
(5) and (77) as 

Ch-^ j \V{v)-f'/^dv<C . 

d(v)>2R 

Together with (88), (89), (102), and (103) this gives the proof of (10), and therefore finishes 
the proof of (7). 

Step 3: Properties of the density. 

We will now show that the density matrix 7 in (101) satisfies the two requirements (8) 
and (9). 

The density of 7 is 

M 

Pli^) = Y] 0r,k{xfpkix) + / ipl{x)4>l{x)pu{x)i{u)-^ du , (104) 

k=i Jd(u)<2R+l 

where pfe for = 1, . . . , M is the density of the density matrix 7/; and pu for u G is the 
density for 7^. We first control the 6/5-norm and the 1-norm we 

use the combined Daubechies-Lieb-Yau inequality (Theorem 16) with a = l3^/'^h~^ < 1/2, 
V = 2zk, and m = h~'^ to obtain that 

> Tr[e,,fe7fce,,fe(r^(-iW) - Zk\x - rk\-^)] > ^Tr[er,klkOr,kT^{-ihV)] 

- Cz^hr^ Ch-\l'\^'^ - Czlh' 
> \ TV {er,kikQrM-'^h^)\ - ch-H-^'^ , 

where the constant C depends on z^. Hence we have that 

Ti[Tp{-\hV)er,klkOr,k] < Ch-^S-^/^ = Ch-^/^ . (105) 
Using (14) with a = < 1 and m = (105) implies that 

(\ 18/25 



(\ 18/20 



^3/10 



< C/i-=^6/25;^-9/5^21/25+ ^^-9/4^3/10 ^c-^-WS^ ^^gg) 

where we have used that r = h and that h is bounded above by a constant. Likewise we find 

9l,pk<Ch-'/\ 

The case when 1/2 < /3^/^h~^ < 1 is more complicated. We have to treat the region 

within the radius r_ = /i^ from the nucleus z^ differently. Let 0±{x) = 0±{\x — r^l/Zi^). 
Using the relativistic IMS formula (Theorem 13) and Theorem 14 with £ = m = h~'^, 

a = ^V2/j-i^ and Q = {\x - rk\ < 3h'^} we find that 

> Tr[er,klkOrMTp{-^hV) - Zk\x - rkt^)] 

> TV[L7fcL(r^(-iW) - Zk\x - rfc|-i - Ch-^)] 

+ Ti[er,ke+lker,kO+{To{-ihV) - Zk\x - rk\-^ - h-\n)] - ch-^ . 
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To treat the first term we use the inequahty (see (21)) 

7r|x| 

which holds for aU < s < 1/2 and As,Bs > being constants depending only on s. Hence, 
using that h is bounded above by a constant and that 1 < /3~^/^/i < 2 we get 

> Tr[e_jk0-{T(s{-ihV) - Zk\x - n]-^ - Ch'^)] > 'IV[^_7J_(A,(-A)^ - C/i'^)] . 

We appeal to the standard (Daubechies)-Lieb-Thirring inequality 

TV[(-A)^^_7fe^_] > c j (^^p,)(=^+2^)/3^ 

which holds for all s G (0, 3). We obtain that (with all constants depending on < s < 1/2) 

Tr[9^^kO-{Tp{-ihV) - Zk\x - r^r^ - Ch'^)] > c J (^^p,){3+2.)/3 _ J (^2 

> (c/2) J (^2 p^)(3+2.)/3 _ . 

Using the Daubechies inequality (Theorem 9) we find as above that 

Tr[er,kO+lkOr,k0+{Tfs{-ihV) - Zk\x - rk\~^ - h'^xn)] 
> cTr[er,k9+lk0r,k9+T/3{-ihV)] - Ch-^/\ 

By choosing s sufficiently close to 1/2 and using that h is bounded by a constant we conclude 
that 

0>cj (^2 p,){3+2s)/3 ^ cTT[er,k0+lker,k0+Tl3{-ihV)] - Ch~^/^ . 

As above it follows from this, choosing s sufficiently close to 1/2, that we still have 

I {el,pkf/' < ch-^'/\ j el,pk < ch-^/' . (107) 

Using that r = h and that from (5) |l^(x)| < Cd{x)~^ we also have 

I {h-X^\v.\^/^f/^ < ch-'''/^ , I h-^el^\v.ff^ < ch-^''' . (108) 

We move to the second term in (104). By the rescaled semi-classics (Theorem 23) we have 
on the support of ^u<P+ that (for /(n), sec (95)) 

p„(x) - 2V2(37r2)-i/i-3|F(a;)_|3/2| < C7/i-2-ViO/(^)2i/io^(^)-9/io + c/i-2|F(x)_|3/2 , 

where we have used that on the support of ipu4>+ we have |F(x)| < Cd{u)~^ < Cr~^ < 
Ch~^ < Chj3~^, since d{u) > r/2 if ipu4'+ is non-vanishing. We moreover have on the 
support of Lpu<j)+ that |y(.x)_|-^/2 < Cf{tif < C f{uf^l^^ iiuy^/^^ . For r/2 < d{u) < 1 this 
is because l{u)~'^ > d{u)^^ = f{u)'^ > f{u) (see (61)) and for d{u) > 1 we simply use that 
i{u) < 1 and f{u) < 1. Hence 

||.^^02+(^^_2l/2(3^2)-l;^-3|^_|3/2)||^^^ < C/.-^-Vl0;(^)21/10^(^)8/5 ^ (lOg) 
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Using (101) and (102), we have that 

||p-,-2^/2^37r2)-i/i-3|y_l3/2|' 



31 



16/5 



(110) 



M 



< 



T.i\KkPkk/,+ch-^di,\v_\'/%/,) 



k=l 
+ 



ld( u 

I 



d{u)<2R+l 



|<^^</,^(p„-2V2(37r^)-^^-=^|F_|3/2)||^/^£(n)-=^d^ 



d{u)>2R+l 



Ch-^^l<t>l\V-f\^^i{u)-^du. 



The last term is non-zero only in the case = in which case it is easily seen by (6) and 
(77) to be bounded by C^-^/^. Thus, combining (106) -(109), (110) implies that 

\\p^ - 2V2(37r2)-i/,-3|y_|3/2|| < ch-^ + C/i-2-Vio f f{uf'/'^eiu)-y^du. 

' Jc-^r<d{u)<2R+l 

The last integral is easily seen to be bounded and we arrive at (9). 
To control the integral of the density we estimate 

J p^{x)dx -2^/\3Tr'^)-^h-^ J \V{x)-\^/^dx 

M 

< T.i\\^lkPk\\i + ch-'\\el,\v.n) 



k=l 
+ 



[ [ ^l4>l{Puix) - 2V2(3vr^)-i/,-3|y(^)_|3/2) 

Jd{u)<2R+l J 



(u) ^du 



+ 



/ Ch~V^cl>l\V^f/\i{ur'du. 

Jd{u)>2R+l 



As before the last term is bounded by Ch For the middle term we again see from the 
rescaled semi-classics (Theorem 23) that 

^id^iip^ix) - 2y\3n'r'h-'\v{xu'/')dx 

< Ch-'^+^'^f{uflH{ufl^ + Ch-^ j ifl{x)(t>l{x)\V{x)-fl'^dx 

where we have used that (3 <h^. The estimate (8) follows since both integrals 

f{ufl^i{ufl^i{u)-^du and / f{ufdu 



/■ 



are bounded (recall that f{u) is given in (95)). This finishes the proof of Theorem 4, except 
for the continuity of the function S from Lemma 25. We will need a lemma to prove this. 
This lemma also gives an alternative characterization of the function S. 

Lemma 27 (Scott-corrected pushed-up Hydrogen). Let S : [0,2/7r] ^ M &e the 
Junction from Lemma 25. Then there exists a constant C > such that, for all a G [0, 2/7r] 
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and K G (0, 1], 



TV [A/-a-2A + a-4 - a'"^ - |^ 1"^ + _ 

- (27r)~3 J [i/ - \v\-^ + k]_ dpdv - S{a) < Ck^/^^ . (Ill) 

Here, as before, y/—a~^A + — = —A/2, when a = 0. 

Proof of Lemma 27. A simple rescaling, changing x — > k~^itx/2, gives 

TrfV-a-^A + a-4 - a^^ - If + k1 = KTrT^-Z^^^^^A + /3-2 - /J-i - ^ + l] , 

where j3 = kc? and h = 2^}^"^ /-k. We have /? < /i^. 

The semi-classical integral may be rewritten in the same fashion, 

(27r)-3 j [\p^ - \v\-^ + k] _dpdv = k {2iTh)-^ J [\p^ ~ ^ + ^] - '^P^'" ' 

Since the potential V{x) = — 1 satisfies the assumptions of Theorem 4 we see that there 
exists a constant C > such that 

TV W-p-^h?A + /3-2 - p-^ -JL. + i\_ 

- (27r/i)-'^ / [ip^ - ^ + ll dpdv-h-^^S{a) <Ch-''+^'^\ 
J '-"^ 7r|f| TT^ 

Using that h = 2k^I'^/'k gives (111). □ 

We can now, using the alternative characterization of the function S in Lemma 27, finish 
the proof of Theorem 4. 
Step 4: Continuity of the function S. 

We recall that 



, /3 = 



It suffices to prove continuity of 

Tr[r„2(-iV) - + = Tr[V-a-2A + a-4-a-2 _ 

at all ao G [0, 2/7r], for any ac G (0,1] fixed. Then continuity of S follows from (111) by 
uniform convergence as k — > 0. 

We first prove the continuity at ao = 0. 

Lst x> = X|p|>A! X< = X|p|<A foi' some A > to be chosen below. Note that (Fi — 
r2)(ri - T2Y > implies FiT* + TaT* < TiT* + TaT*. Using this with Ti = ei/2^^|x|-V2^ 
T2 = £^^^'^X>\^\^^^'^ for some e > which we choose later, we get the operator inequality 

ra2(j5) - Ix]-^ + K (113) 

> X>{Ta^{p) - (1 + e'^)\x\-^ + '^)X> + X<(7;,2(p) - (1 + e)\x\-^ + k)x< . 

Here and in the sequel we write T„2(p) for the operator T„2(— iV) (and similarly for other 
operators). Since T^2 > T^2 for ai < 02, and T^a{p) > a^^\p\ — a^^, (113) implies that, if 
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a G (0, ^] for some A > 0, then for all e > 

T^-^{S>)-\x\~^ > x>iA-^\p\-A-^-il + e-^)\x\-^ + K)x> 

+ x<{T^^ip) - (1 + + '«)X< • (114) 

Since \p\ — 2/(7r|x|) > 0, we have that 

l^-i|p|_(l + £-i)|^|-i>0, 
if A < e/(27r), and now assuming £ < 1. Furthermore, for A > 2A''^ we have that 

x>{Ia-'\p\-a-^)x>>o- 

This implies that, if e < 1, A > 2A-^, a e (0, A] and A < £/(27r), then by (114) 

T^2{p)-\x\-^ + K > x<(ra2(p)-(l+£)|x|-^ + /c)x<. (115) 

Since, by a Taylor-expansion, T^2{p) > Tq{p) — (ap^)^/8, and since x< = X|p|<A) we have 
that, stih for a G (0, A], 

r,2(p) - \x\-^ + K> x<{To{p) - a^AVs - (1 + e)\x\-' + k)x< . (116) 

Let 

1a,K = X(-oo,0] (Ta^iP) - + «) • 

Then (116) and the fact that Tq > T^-z imply that, for a € (0, A], 

Tr[ro(p) - \x\-^ + > TV[T^2(p) - \x\-^ + k]_ = Tr [7„,4r„2 (p) - \x\-^ + k)] 

> Tr[7„,«x<(7o(p) -a2A78-(l + 6)1x1-1 + /t)x<] . ^^^^^ 

If K G (0, 1], a G (0,^], A > 2A~^, and A < l/(27r) we will show the a priori estimate 

1V[7a,«X<] < and TV[7„,«x< l*r'x<] < C^-^/^ (118) 

The combined Daubechies-Lieb-Yau inequality (32) gives that for positive constants Ci,C2 
such that a < 2/(Ci7r), we have 



Tr[T^2{p)-Ci\x\-^ + C2K]_ > -Ca^/^-C [ [Ix]-'^ + k)' 

- Ca^ [ {\x\-^ + Kj^dx > - Ck-^/^ . 



la<\x\<CK- 

If Q G iO,A] and A < l/(27r) then a < 4/(57r) and hence we obtain from (115) with e = 1 
that 

> Tr[T^2{p) - Ixpi + k]_= Tr[7«,4T«2(p) - \x\-^ + k)] 

> 1V[7„,^x<(7;2(p)-2|x|-i + «)x<] 

5/2 11 K 

> - Ck-^/^ + ^Tr[^a,nX<m-^X<] + f 1V[7a,«X<]- 
This gives (118). 
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Choose A = 2A-\ A = e/(27r). We combine (117) and (118) and use the variational 
principle to conclude that for a G (0, e/(27r)], e < 1, and k G (0, 1], 

Ti[To{p)~\x\-^ + n]_ > Tr[r„2(p) - |xr^ + 

> 'IV[(x<7a,«x<)(7b(j5) - \x\-^ + K)]-CK-^l\a^e-^ + e) 

> Tr[ro(p) - \x\-^ + - C^-^/^^aV^ + e) . 

Finally choose e = a^/^; then a < {2tt)~^^^ implies that a G (0, e/(27r)] and e < 1. 
Therefore we have proved that for any a < (27r)-5/3 and k G (0,1], 

Tr[To(p) - + > Tr[T^2(p) - Ixj-i + 

> TV[To(p)- + -C«;-3/V/^ 

which proves continuity from the right of Tr [Tf^2 (p) — \x\~^ + k] _ at ao = for any k G (0, 1] 
fixed. (Notice that the above has not been optimized in k.) 

We now prove the continuity at any ckq G (0, 2/7r). Note first that, for < ai < Q2, 

Te^iip) > T^2{p) > {a^'aifT^2{p) . (119) 

Assume first that a > ao, and let 'fa^K be defined as above. Then, using (119) and the 

variational principle, 

TV[T^2(p) - \x\-^ + k]_ < Tr[T^2(j5) - \x\-^ + k]_< Tr [7,,4r„2 (p) - {x^' + k)] 

< Tr[7„,«(r„2(p) - + + [{aa^'f - 1] Tr [7„,«T,2 (p)] 
= TV[r<,2(^) - + k]_ + [(aao ')' - 1] TV [7„,«r„2 (p)] . 

It remains to show that [(aag "^)^ — 1] Tr [7q^kT'q,2 (p)] ^ as a — ^ ccq. For this, it obviously 
suffices to show that Tr [70, ^2^02 (p)] is uniformly bounded for, say, a G (ao,A] for some 
A G {ao,2/7r). But this follows as in the proof of (118). This proves continuity from the 
right of Tr[Tg,2{p) — \x\~^ + at ckq € (0, 2/7r). To prove continuity from the left, assume 
a < ao, and let 7ao,K be defined as above. Then, by (119) and the variational principle, 

Tv[r,2(p) - \x\-^ + k]_ > Tv[r,2(p) - \x\-^ + k]_ = iv[7,„,«(r,2(p) - \x\-^ + «)] 

> TV[r„2(p) - |:r|-^ + + [1 - (ao«"')'] IV baoA^ag (p)] • 

As before, the last trace is finite by arguments as in the proof of (118) (since oq < 2/7r). 
This proves continuity from the left, and therefore, continuity, of Tr[rQ,2(p) — l^l"^ + 
at ao G (0, 2/7r). 

Finally we prove the continuity at ao = 2/7r. Here, arguments as in the proof of (118) 
are no longer at our disposal. Therefore, let e > 0, and let ^ao,K be defined as above, and 
choose ^1, . . . , € Cq°(R^), {(f)i, = Sij, such that 

Tr[jM{T^2{p)-\x\-' + K)] (120) 
< Tr[7„o,4T,2(p) - |x|-i + k)] + e/2 = Tr [r,2(p) - \x\-' + n]_ + e/2, 

for jNixjU) = J2jLi4'j{x)4'j{y)- This is possible since the operator is defined as the 
Friedrichs extension from Co°(M^). (Here, both and the ^j's depend, of course, on e). 
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Recall that is finite dimensional and (pj e C^(M ). Using this, (119), and the variational 
principle gives that (for any a G (q;o/2, ao)), 

T^[7iv(T,§(p) - + Ac)] = Tr[^NT,2{p)] + 1V[7iv( - l^l"^ + k)] 

> Ti[jM{T^2ip) - \x\-^ + k)] + [(ao - 1] T^ilNTM] 

> Tr[T^2{p) - \x\-^ + k]_ + [(ao ^a)2 - 1] TV[7^r,2/2(p)] . (121) 
Choose now S > such that 

ae (ao-5,ao)n(ao/2,ao) [{aQ^af - 1]Tt[jnT^2/2{p)] > - e/2 . (122) 

Then, combining (120), (121), and (122) (and using (119) again) we have proved that, for 
all e > there exists < S < ao/2 such that 

a G (ao - S, ao) 

IV[r„2(p) - \x\-^ + i^]_ > Tr[T^2{p) - + Ac]_ > IV[r„2(p) - + - e. 

This proves the continuity from the left of Tr [rQ,2 (p) — at ao = 2/7r, and therefore 
finishes the proof that S : [0, 2/7r] — s- R is continuous. 

This completes the proof of Theorem 4. □ 

5. Local relativistic semi-classical estimates using new coherent states 

In this section wc study the sum and the density of the negative eigenvalues of the localised 
Hamiltonian (pHgcj), with compactly supported and Hjs = T^(— i/iV) + V{x). Here, is 
given by (3), and F is a (sufficiently) regular potential (see below for details). For the most 
part we suppress the index P but all estimates, in particular the constants C, will be uniform 
in/3G [0,1]. 

We first recall the definition and the main properties of the coherent states (operators) 
introduced in [36], where all proofs can be found. These coherent states are denoted by Qu,q- 
Let 1/a > ^ > 0. The kernel of Qu,q is given by 

Gujx, y) = (vr/i)-"/2e-'^(^-")'+^«(^-^)/'^-3^(^-^)' . (123) 
A first important property of these operators is their completeness. 

Lemma 28 (Completeness of new coherent states). The coherent operators Qu,q satisfy 

where x denotes the operator multiplication by the position variable x. Here Gb{v) = 
(6/7r)"'/^e~^'' with b = 2a/ {1 + /i^a^). Note that Gb has integral 1 and hence 

We shall consider operators of the form 

Gu,q f{Au,q) Gu,q dudq , (126) 



I' 



where / : R — >^ R is any polynomially bounded real function. As we shall see in the next 
theorem the integrand above is a trace class operator for each {u,q). The integral above is 
to be understood in the weak sense, i.e., as a quadratic form. We shall consider situations 
where the integral defines bounded or unbounded operators. 
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Theorem 29 (Trace identity). Let / : M ^ M and V iR"' be polynomially bounded, 
real-valued measurable functions and let 

A = Bq + Bix - \hB2V 

be a first order self-adjoint differential operator^ with Bq G M, .81^2 € M". Then 
Gu,q f{-^) Gu,qV{x) is a trace class operator (when extended from Cq^{W")) and 

Tr[gu,gf{A)g^,gV{x)] = J fiBo + Biv-hB2p)Gb{u-v)Gb{q-p)G^hH)-^iz) 

X V{v + h^ab{u — v) + z) dvdpdz . 

In particular, Tr[C/,^ = 1. 

We shall also need the following extension of this theorem, where we however only give 
an estimate on the trace. 

Theorem 30 (Trace estimates) . Let /, A be as in the previous theorem. Let moreover 

(f) G C"'~''^(M"') be a bounded, real function with all derivatives up to order n + 4 bounded, 
and let V,F ^ C^(M") be real functions with bounded second derivatives. Then, for h < 1, 
1 < a < 1/h and b = 2a/(l + /i^a^) we have, with a{u,q) = F{q) + V{u), that^ 

Tr [g^^q f{A) 0(x) (F(-iW) + V{x)) </.(x)] 
= j fiBo + Biv + B2P) Gb{u - v)Gbiq - p) 

X {(l){v + h^ab{u — v)Y + -Ei(n, v)) a{v + h^ab{u — v),p + h^ab{q - p)) 

-\- E2{u,v;q,p) dvdp , 

with ||£^i||oo, ||-£'2||oo ^ Ch%, where C depends only on 

sup ||9'^0||oo, sup ||9^F||oo , and sup ||5'^F||oo • 

\v\<n+i \v\=2 \v\=1 

(Note that the assumption 1 < a < 1/h implies 1 < b < 1/h.) 

We will use the above theorem to prove an upper bound on the sum of eigenvalues of the 

operator F{—ih'V) + V{x), in the case when F{q) = Tp{q) from (3) with (3 £ [0,1] (equal 
to ^JlPW^nf^ - (3^^ for (3 G (0, 1], and to \q^ when /3 = 0). This is done in Lemma 34 
below by constructing a trial density matrix on the form (126). 

To prove a lower bound on the sum of the negative eigenvalues one approximates the 
Hamiltonian F(— iW) + V{x) by an operator also represented on the form (126). 

Theorem 31 (Coherent states representation). Consider functions F,V ^ C^iW^), 
for which all second and third derivatives are bounded. Let a{u,q) = F{q) + V{u), then for 
a < 1/h andb = 2a/(l + /i^a^) we have the representation (as quadratic forms on Cq°{M."-)), 

F(-iW) + V{x) = I g^^qH^^qg^,^-^ + E , 

with the operator-valued symbol 

Hu,q = (t{u, q) + 35Acr(n, q) + 5„cr(u, q){x - u) -\- dqa{u, q){-ihV - q) . (127) 



•^The operator A is essentially self-adjoint on Schwartz functions on R". 

"^The operator f{A) Gu^q d}{x) (F(-i/iV) + V{x)) 4){x) is originally defined on, say, Co°(R"), but it is 
part of the claim of the theorem that it extends to a trace class operator on all of L'^(M"). 
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The error term, E, is a hounded operator with 

||E|| < W<^\U + ChH \\d"(^\\oo ■ 

\u\=3 \u\=2 

Let us recall our convention that a;_ = minja;, 0} and that x = X(-oo,o] denotes the 
characteristic function of (— oo,0]. 

The next theorem is the main result of this section. 

Theorem 32 (Local relativistic semi-classics). For n > 3, let CP e Co"+^(M") be 
supported in a ball Bi C M" of radius 1 and let V G C^(i?i) be a real function. Let 
< (3 < I, h > 0, and let afs{u,q) = T(s{q) + V{u) and H/3 = r^(-i/iV) + V{x) with 
Tpiq) = VrVT^ - for P G (0, 1] and To(g) = ig^. 
Then 

T:t[<J)Hp(P]_- {2TThy j (P{ufap{u,q)-dudq < 0/1""+^/^ 

The constant C > here depends only on \\(j)\\cn+A, \\V\\cz,^ and the dimension n, but not 
on (3£ [0,1]. 

The important property for our method to work is that the second and third order 
derivatives of the kinetic energy function Tp{q) are bounded uniformly in q and (5. Thus 
the error term above is independent of /3 G [0,1], and in particular the same as for the 
non-relativistic case, — h'^A/2 + V, which corresponds to the limit /3 — > 0. We prove upper 
and lower bounds and start with the lower bound. 

Lemma 33 (Lower bound on Tr[0if^^]_). Under the same conditions as in Theorem 
32, 

Tr:[(pHp(t)\^ > (27r/i)-" j (t){ufap{u, q)- dudq - Ch'^'+^l^ . 

The constant C > here depends only on ||^||cn+4, ||l^||c3; the dimension n, but not 
on (3e [0,1]. 

Proof. Since (p has support in the ball Bi we may assume without loss of generality that 
V G Cq(M") with the support in a ball B2 of radius 2 and that the norm HFUca refers to 
the supremum over all of M". Wc shall not explicitly follow how the error terms depend 
on ||(/>||c7n+4 and ||T^||c73. All constants denoted by C depend on ||(/>||(^n+4, ||F||(73, and the 
dimension n but, in particular, not on /3. 

We use the Daubechies inequality (Theorem 9) to control various error estimates. Since 
Tf^iq) > Ti{q) for /? G [0, 1] we may use it with /3 = 1. Then, uniformly in /3 G [0, 1], 

Tv[cl>Hf,cl>]_ > I a,{u,q).^^> -Ch--. 

ueBi 

Consider some fixed < r < 1 (independent of h and (3). If h> t then we get that 



dudq 



Thy 

We are therefore left with considering h < t. 



^We use the convention that = sup|j,|<j, HS^i/'llc 
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If we now use the inequality [x + ?/]_ > + [y]_, which we will do frequently without 
further mentioning, and Theorem 31, we have that 

+ TV [(t){e^-p-^h^A + /3-2 - e/^-^ - (7(6-^/^ + h^b))(t>] _ . (128) 

Here, < £ < 1/2 and 

Hi^^q = ^(.u, q) + j^Aa{u, q) + q){x - n) + dqa{u, q){-ihV - q) 

with a{u,q) = (1 — £)Tg{q) + V{u). The second trace can be estimated from below using 
the Daubechies inequality (Theorem 9) with a = P^^^h~^, m = h~^. Then 

Tr[<l){ey/-p-^h^A + p-^ - e^^ - C{b-'^/^ + hH))(f>]_ 

= £TV[(/)(V-/?-^/i2A + /3-2 - p-^ - Ce-\b-^/^ + h^b))(l)]_ 

> - Ceh-'' I (£-1 (6-3/2 ^ ^2^)) l+n/2 



(129) 



We shall eventually choose e =\(b 3/2 _)_ /j25^ Note that then e < 1/2, and that the bound 
in (129) is - Ch-''{b~^/'^ + h%), uniformly for /3 € [0, 1]. 

By bringing the negative part inside in (128) we obtain the lower bound, 

We first consider the integral over u outside the ball B2 of radius 2, where = 0. Using 
Theorem 29 (with f(t) = and V replaced by cp'^) and j (fP' < C, we get that this part 
of the integral is 



{l-e) 



Tp{q) + (n + (n - l)/3(?2)/[46(l + /3g2)3/2] + ^ . _ g)/yTT^ 



dudq 



xGb{q- p)Gb{u - v)G^hH)-^ iz)4>iv + h ab{u - v) + z) dvdpdz ^^^^^ 
> (1 — e) J (t'iz)^ J G]){u — v)G(^f^2i,yi{v + h?ab{u — v) — z) dudvdz 



U^B2 



dqdp 



(1 — e) J (j){z)^ J Gi,{u)G^fi2f,^-i{v — z) dudvdz 



zeBi 



{l-h?ah)u^B2-v 



[Tp{q)+q-p/Vl + P^]_G,{p) 



dqdp 
(27r/i)" 



The integration over u, v is obviously bounded by 1. In fact, the ^-integration can be shown 
to be exponentially small, i.e., less than C e~^^ , but this will not be necessary. 
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The domain of integration for the variables q,p is contained in the set {{q,p) \ \q\ < 2|]9|}. 
Then, 



J [Tp{q) + q-p/^/l + ^]_G,{p)dqdp 



> -C 



I MH Q^^^ ^^^^ >-C I br+2G6(p) dp= -C 6-("+2)/2 _ 



l?l<2bl 



It follows that the integral over n S2 is bounded from below by — C^~"6~^/^, since h> 1. 

For the integral over « G S2 we use Theorem 29 as before. This time, expanding to 
second order in z at the point 2; = and using the crucial fact (which we shall use without 
mentioning later) that, for any A > 0, 



j Xj Gx{x) dx = 0, J Gx{x) dx = C X'""/^ , 

implies that 

Tr[(PHf}(j)]_ > J [<p{v + h^ab{u - v)f + ChH] Gb{u - v)Gb{q - p) 

dudq 



(130) 



U&B2 



(27r/i)' 



;dvdp - Ch-'^ih-^/^ + hH) , (131) 



where 



Hi%v,p) = a{u, q) + ^A(t(u, q) + dua{u, q){v - u) + dqa{u, q){p - q) . 



The rest of the proof is simply an estimate of the integral in (131). This analysis is an 
elementary but tedious exercise in calculus. For the convenience of the reader it is given in 
detail in Appendix B below. □ 

Lemma 34 (Construction of a trial density matrix). Under the same conditions as 

in Theorem 32 there exists a density matrix 7 on L^(R"') such that 



Tr[(f>{Tp{-ihV) + V{x))4>j]< J Hufapiu,, 

Moreover, the density of 7 satisfies 



dudq 
(27r/i)" 



p^{x) - (27r/i)-'^w„|y_|'*/2^2 + /3|y_|)"/2(x) < c/i-'^+^/io , 
for (almost) all x E Bi and 



(132) 



(133) 



xfp^{x)dx-{27rh)-''u;n / <p{xf \V-\''/\2 + P\V-\r/\x) dx 



< C/i-"+6/^ , (134) 

where ujn is the volume of the unit hall Bi in R" . The constants C > in the above estimates 
depend only on n, \\4>\\Qn+4, and \\V\\c3, but not on P E [0, 1]. 



It is convenient to introduce the function 

fOO 



/•oo 

v{t) = n / xiTpip) + 1] \pr' d\p\ = |t_r/2(2 + p\t^\r 

Jo 

(Recall that x is the characteristic function of M_.) 



/2 



(135) 
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Proof. We will occasionally drop the index P in if^ and a p. It is important to realize, 
however, that all estimates are uniform in /3. We first note that since Ti{p) < Tp{p) < 
To{p) = we have that 

\p\-C<a{v,p)<^p'' + C. (136) 

Let us start by choosing some fixed < r < 1. For h > t and for some C > we have by 
(136) that 

/^MMn,.)-^ + Cr-/^^— «/^>0, 
and that for any s > 0, 

(27r/i)-"?7(1/(x)) < Cr-^/i-"+" . 

If h > T we may therefore use 7 = 0, and s = 9/10 and s = 6/5 for (133) and (134), 
respectively. Prom now on wc assume that /;, < r and, if necessary, that r is small enough 
depending only on (/) and V. Also, as for the lower bound, we may assume that V € Co(M"') 
with support in the ball -B3/2 concentric with Bi and of radius 3/2. 

In analogy to the previous proof for the lower bound we define now for each [u, q) an 
operator hu^q by 

a{u, q) + -^Aa^u, q) + Va{u, q) ■ {x — u, —ihV — q) if n G -B2 

if u ^ ^2 ■ 



hu,q 



hu,q{v,p) 



The corresponding function is 

a{u, q) + ^Au{u, q) + Vcr(n, q) ■ {v — u,p — q) if u e B2 

if u ^ S2 ■ 

As for the lower bound we shall choose a = h~^/^; then a < h~^. In fact, we will assume 
that (1 — h'^ab) > 1/2. Recall here that b = 2a/(l + /i^a^) (i.e., in particular a <b < 2a). 
Similar to (172) (for e = 0) we have for u e B2 that 

\hu,q{^^P) -<^{v,p) -iv,p{u-v,q-p)\ 

< C\u - v\{b-^ + \u- v\^) + C\q - p\{b-^ + \q- p\^) , (137) 

where 

^v,p{'^, q) = iE^cr{v,p) -\'Y^ didjTfi{p)qiqj - ^ ^1 9idjV{v)uiUj . 

Recalling that x is the characteristic function of M_ we define 

1 = 1 Gu,q X [hu,g\ Gu,q ■ (138) 

Since < x[K,q] < 1 it follows from (125) that < 7 < 1. 

We now calculate Tr[70iJ^^] = Tr[70(r^(-i/iV) + V{x))(p]. From Theorem 30 we have 
that 

Tv[70(r^(-iW) + y(x))(/>] 

= J x[hu,q{v,p)]Gbiu-v)Gb{q-p) E2{u,v;q,p)+ (139) 
((j){v + h?ab{u — v))"^ + Ei{u,v))a{v + h^ab{u — v),p + h^ab{q — p)) f^f^ dvdp , 
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where Ei,E2 are functions such that ||-Ei||oo + II-E2II00 ^ Ch^b. The rest of the proof of 
(132) is a tedious, but elementary analysis of this integral. A detailed analysis is presented 
in Appendix B below. 

It remains to estimate the density together with J (p{x)'^p^{x) dx. By Theorem 29 and 
(138), 7 is easily seen to be a trace class operator with density 

x[K+v,q+p{v,p)\Gb{u)Gb{q)G{h2^)-i{x -v- ahu) dvdp -^^-^ . (140) 

The proof of (133) and (134) again relies on a detailed analysis of this integral. As for 
the estimate on the energy above this analysis is an exercise in calculus. Although it is 
still elementary this analysis is more complicated than in the case of the energy. For the 
convenience of the reader the analysis is given in detail in Appendix B below. □ 

Appendix A. Various Proofs 

In this appendix we collect proofs of various results mentioned in Section 2. 

Proof of Theorem 11 (Operator inequality critical Hydrogen). Let / G 5(M^) and t > 
(to be chosen below). By Schwarz' inequality. 



^ i.r f li^ifitK,^. ,,41) 

We first compute the integral in q. Since {\q\^ + l^l*)"^ < \q\~'^ — |g|*~^ + we get 

Jm:^ \p - q\ \q\ + \q\ JRs\p-q\^ 

Note [19, 5.10 (3)] that, for < r,c7 < n, with < t + a < n, 



I 



\y - ^r-^i^r-^^ dz = ^^-^-^^^g'^ iyr+^-" , (142) 

where = 7r~'^/^r(T/2). In particular, if n = 3, then 

/ \y-z\-^\z\-''dz = kr\y\^-'^ for re (1,3), (143) 

with 

'~ r(^)r(i) • ^''^^ 

It follows that, for 3 < 2t < 5, 



/ 



1 \p\^ + \p\\^^ ^ k2\p\ + {k2-h-t)\pt' (145) 



/r3 |p- IqP + 

+ {ke-2t-h-t) \p\''-' + ke-2t \pf'-'. 

We see from (144) that k is symmetric with respect to r = 2. Using T{1 + z) = zr{z) in 
the denominator in (144) with z = 1 — r/2 and the relation r(z)r(l — z) = tt/ sin(7rz) (for 
< z < 1) in the denominator and numerator we obtain 

2 tan(7rr/2) 
= l-r/2 ' 
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which shows that k is decreasing on (1,2) and increasing on (2,3). 

Hence from (145), choosing further t > 5/3, we find, for positive constants ^(t_i)/2, 
B{t-i)/2, that 



/ 



J2 I Ult 



^ dq < k2\p\-7r^A^t-i)/2\pt'+7^''B^t-i)/2- (146) 



Since k2 = vr'^, this and (141) imphes that 



-/ ^-^^dx < f \M{\p\-At-i)Mpt' + Bit-i)/2)dp, (147) 

which imphes the operator inequahty, for aU t G (5/3,2), 

V^-^ > ^(t-i)/2(-A)(*-^)/^-S(t_i)/2. (148) 

Choosing t = 2s + 1 proves (21) for s e (1/3,1/2). For s G [0,1/3], (21) follows from the 
existence of positive constants ^(r-i)/2) -S(r-i)/25 given r G [1, 5/3], t G (5/3, 2) and positive 
constants ^(t-i)/2) -B(4_i)/2, such that 

^(t-l)/2W*"^--B(t-l)/2 > ^(r-l)/2br"^-^(r-l)/2- 

□ 

Integral representation for the relativistic kinetic energy. We shall here give a 
self-contained presentation of the integral formulas for the relativistic kinetic energy. The 
relativistc kinetic energy will be given in terms of the modified Bessel functions of the second 
kind, K,y. To identify the modified Bessel functions we use that [1, 9.6.23] 

roc p—wt 

Mt)= —==dw, t>0, (149) 
Jl Vli^ — 1 

and the recursion relation [1, 9.6.28] 

K,+i{t) = -t''^(t-''K,{t)), t>0. (150) 
dt 

We emphasise that we use these properties only as definitions of the Bessel functions, and 
derive all other properties of these functions that we need. Note that Ki, : R_|_ — M are 
smooth functions. 

Consider the function G L^(]R") (the Yukawa potential) whose Fourier transform is 

G-(0 = (27r)-'^/2^|C|2 + m2)-i. 

Using that v^^ = e""'' du we get from the Fourier transform of Gaussian functions the 
following integral representation for G, 

POO 

G;^(z)= / (47ru)-"/2e-"^'«-I^P/(4n)^^_ (151) 

^0 

It follows from this that G is non-negative, smooth for z ^0, and indeed in L^(M"). 

For odd n the above integral can be explicitly calculated. For even n it is as wc shall now 
see expressible as a modified Bessel function Ki, of integer order v. By a simple change of 
variables {2w = v + with v = 2mu/\z\) in the integral (149) we see from (151) that 
G^{z) = {2'k)~^ KQ{m\z\) . From the recursion formula (150) we then find inductively that 
for even n 

G^{z) = m("-2)/2(2^)-"/2|^|-('^-2)/2^^^_^^^^(^|_^|) _ (152) 
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In fact, the same formula holds for all n, but we do not wish to discuss the modified Bessel 
functions of fractional order (one could simply take this formula as their definition) . 

Lemma 35. The heat kernel for the operator V— A + on L^(M"') is given by 

exp(— A + m'^){x, y) 
= -20^^- i(x-2/,t) 

for t > 0. 

Proof. It suffices to show that the two tempered distributions on M"+^, 
^exp(-|t|V-A + m2)(x,0) and - 2dtG^+^{x,t) , 

have the same Fourier transform. The Fourier transform as a function of ^ = {p, s) with 
p G M" and s G M of the first distribution is 

J — 00 Jo |p| + S + 777. 

The Fourier transform of the second distribution above is 

-2is(2^)-(«+i)/2 



-2isG-i(p,s) = 



The last identity in the lemma follows from (150) and (152). □ 

If we set a; = y in the above lemma we find the following integral formula for the modified 
Bessel function 



1 / t 



(n-l)/2 



For n = 3 this simplifies to 

K2{t) = t r e-*^^s2 ds (153) 
Jo 

from which we immediately get the estimate 

K2{t) < Ct-h-^/^ . (154) 

Proo/ 0/ T/ieorem i5 (Relativistic IMS formula). By scaling, it suffices to prove the 
statement for a = 1. We start from the identity 

(/,(V-A + m2-m)/) = J \f{x)- f{yfF{x-y)dxdy 

with 

F(.-,) = ^ ^^H-7l) , (155) 
47r"^ \x — 
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where K2 is the modified Bessel function of second order defined above (see (149)-(150)). 
The identity follows from Lemma 35 (for a proof, see [19, 7.12]). Then, 



(/,(V-A + m2-m)/) 
= / \f{x)-f{y)\^F{x-y)dxdy 

^ [eu{xf\f{x)\^ + eu{yf\f{y)\^] F{x - y) dii{u)dxdy 

M 

I I [- \{eu{xf + Ouiyf) + 0uix)duiy) - 0u{x)eu{y)] 

J JM 



X [ f{x)f{y) + f{x)f{y) ] F(x - y) diJ,{u)dxdy 

= I f \eu{x)fix)-euiy)f{y)\^F{x-y)dii{x)dxdy 
J JM 

+ [ [ [-kiGuixf + 0uiyf) + 0uix)euiy)] 

J JM 

X [Jix)f{y) + f{x)J{y)]F{x-y)dn{u)dxdy 
= [ (^„/,(^^A+^-m)^„/)d/x(t^) 

JM 

- [ [ {eu{x)-eu{y)Yf{x)J{y)F{x-y)dii{u)dxdy. 
J JM 

This proves (23) with L given by (24)-(25). We now show that \\Lq^\\ < Cm~^||V6'„||^ for 
fixed u. By (25), Young's inequality, and (154), 

777^ f 

\{f,LeJ)\ < -^WVe^Wl J \f{x)\\f{y)\K2{m\x-y\)dxdy 

/•OD 

< Cm^\\V9u\\lo\\f\\l / t^K2{mt)dt 
Jo 

= Cm-i||V^„||^||/||i. (156) 
This proves that Lg^ is a bounded operator. □ 

Proof of Theorem I4 (Localisation error). Again, by scaling, it suffices to prove the state- 
ment for a = 1. With XQ. the characteristic function of (and L = Lq) we have from the 
representation (25) of L, since is constant on OF, that 

L = xnLxn + (1 - Xn)Lxn + XnL{l - Xn)- (157) 

If Li, T2 are bounded operators, then (Li - r2)(ri - T2)* > implies that TiTl + T2TI < 
TiTl + T2TI. Using this with Li = e^/'^xn,'i^2 = £~^/^(l - Xn)L for some e > which we 
choose later, we get 

L < xnLxn + exn + £"^1 - Xn)LH^ - Xn) ■ (158) 
To bound the first term on the right side recall that ||L|| < Cm~^ ||V^||^ (see (156)). 
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Let us now look at the third term in (158). Since 9 is constant on Q'^ and 
dist(r2'^, supp V0) > using (25) gives 

Tr[(l-xa)i^'(l-Xn)] = / L{x,y f dxdy 

Jxeii'',yen,\x-y\>e 

< Cm^llV^II^ [ K2{m\x - y\f dxdy . 

J xeCl'^ ,y&Q,,\x-y\>i 

Using (154), 



/ K2{m\x — y\)'^ dxdy 

Jxe^'',ym,\x-y\>t 

< Ce"™^ f {m\x - j/|)-^ dxdy 

f oo 



poo 

= Ce-"'^\n\ / {mt)-H'^dt = Cm-^ime)-h-'^^\n\. 



This gives the bound 

Tr[{l-xn)L\l-xn)] < Cr'\\Ve\\t,e-^'\n\ . 

Finally, we choose e = m^^||V0||^. Then by the above the two first terms in (158) are 
bounded by Cm~^ ||V^||^xr2, and the trace of the third term (which we denote Qg) is 
bounded by Cm^-V"^^ II V^||^|0|. □ 

For the proof of the combined Daubechies-Lieb-Yau inequality (Theorem 16) we need the 
following inequality [5]. 

Lemma 36. For f eS{R^), 

/ ^ |/(x)|2 dx < 5^=—- (/, ( V-A + m2 - m)f) . (159) 

Proof. Let = m^vr^^. Then 

f e^^l^l^ 2 ^ f f ~ 1 

^= / I I l/(a^)l dx = —^ / f{pi)- 12 9{P2) dpidp2 , 

m 2tt^ J^3 J^3 \P1-P2r 

with g{x) = /(x)e^^l^'l^. Writing g{p2) explicitly as the convolution with the Fourier trans- 
form of e~''l^l and then applying the Schwarz inequality we get that 

1 f f f |/(pi)pe-l^lV(^^) \Pi\\ ^jp^^^_ 

~ 167r2(7r/i)3/2 7^3 As jRa 1^2-^31-^1^ |^32p 



Since [19, 5.10 (3)] 



/ TZ 1 n^TZ-T^dp2 



\P2 -pi- q\'^ b2p \pi - q\ ' 
we have 

|/(^^)|2e-klV(4.)^^ 



^ - i« 1/2 3/2 / / \ \ \P^\ dpidq. 
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By Newton's theorem [19, 9.7 (5)], 



dq = e-\i\I^^^Uq+ / dq 



\Pi-<l\ \P^\ J\<i\<\pA J\q\>\pi\ 



Jo \Pi\ 



\Pi\ 

Substituting /x = m^Tr~^ we find that 

I < \f{pi)\'^mm{\pi\^,m\pi\}dpi, 

from which the claim follows since VPTT- 1 > (V2- l)min{i2,t} for t > 0. □ 

Proof of Theorem 16 (Combined Daubechies-Lieb- Yau) . We may assume that W{x) < 
otherwise we simply replace W by 

Assume first that va < 3/(167rM). By the Daubechies inequality (17), 

> \ TV[V-a-2A + m2a-4 - ma'"^ + 2t^X{daW<am-i}] _ (160) 
-Cm^l'^ I \W{x)fl'^ dx-Co? i \W{x)\Ux. 

The assumption on the positions of the nuclei implies that X{d^{x)<am-'^} — 
X{\x-Rj\<am-^}-i ^iid so, using the assumption on W , we obtain 

Tr [ V-a~^A + m^a-^ - ma"^ + 2Wx{d^(x)<am-^}] _ 



= Tr[\/-a-2A + m^a-^ - ma~^ - (2z/M|xr^ + CuMma~^)x{\x\<am-^}\_ ■ (161) 
The last equality follows from the translation invariance of — A. By scaling, 

Tr[V-Q-2A + rn^a-^ - ma~^ - {2vM\x\'^ + CuMma~^)x{\x\<am-^}] _ 



TV[V-A + m2 - m - (7|x|-i + C7m)x{|a.|<m-i}] _ , (162) 
with 7 = 2Mi/a < 3/(87r). Using Lemma 36 and the Daubechies inequality, we get that 
Tr[V-A + m2 - m- {7\x\~^ + Cjm)x{\x\<m-^}]_ 



1 — 771 TT 

> (1 - f 7) TV[V-A + m2 - m - 7(1 - f 7)"'( ^^pq + Cm)x{H<^-i}] 



2^-lU|2 



/* /"I — TTl TT uC /I 

> _ 6-7^/2^3/2 / ^ ^^5/2 _ ^^4 / / + ^) dx, 



1 

'|a;|<m~i v/|a;|<m~ 

where we have used that - 1 > 3/8 and 7 < 3/(87r). 
Note that 

/ {\x\~^ + mf^"^ dx < Cm-'^l'^ , j ( ~ + m) dx < Cm , 

J\x\<m-^ -'|a;|<m-i ^ 1^1 ^ 
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and so 



Tr[v'-A + m2-m - (7|x|-i + C7m)x{|,,|<™-i}]_ > -^(t'/^^tV- (163) 
Combining (162) and (163), and using 7 = 2Afz^a, ua < we get that 

Tr[\/-a"^A + m?a-'^ - ma~^ - (2z^M|x|"^ + CuMma~^)x{\x\<am-^}] ^ 

> ~Ca-^{-f^/^ +j^)m> -Cv'l'^a^l'^m. (164) 

Combining (160), (161), and (164), yields (32) for va < 3/(167rM). 
Assume now that va G [3/(167rM), 2/7r]. 

Let e C^{M.^) satisfy < e{x) < 1, 9{x) = 1 for |x| < am-^/4, e{x) = for 
\x\ > am-^/2, (1 - 02)1/2 e C\R), and 

llV^lloo < Ca-V, ||V(l-0^)^/^||oo < Ca~'^m. 

Let Ojix) = e{x - Rj), j = 1,...,M, and OM+iix) = (1 - Yljii ^|)^^^ (the latter is well- 
defined due to the assumption min^^^ — Ri\ > 2am~^). The relativistic IMS formula 
and the localisation estimate (26) used for Qj, j = 1, . . . ,M, being the balls centered at 
Rj with radii 3am~-'^/4 and £ = am~^/4, and f^M+i being the (disjoint) union of the same 
balls and £ = am^^/S, gives the operator inequality 

M 

> 9m+i {V-a-^A + m2a-4 - ma'^ - Cma-^ xn,- + H^(x))0m+i (165) 

M M+l 

+ ^ 0j(V-a-2A + m2a-4 - ma-^ - Cma^^ + W(x)) 0j - ^ Q j , 

with 

Tr[Qj] < Cma-^ 

Here wc have used that OiXn-Oi = 6ijei^, i,j G {1, ...,M}, Om+iX^m+i^m+i = 0, and 

eanM+^o^ < el i ^ M + 1. 

Using the Daubechies inequality on the first term in (165) and the assumption on W in the 
second (noticing that Oj{x)/dii{x) = 9j{x)/\x—Rj\ due to the assumption minj^^^ > 
2am^^), wc get from this that 



Tr[A/-a-2A + m2a-4 - ma"^ + W{x)]- 
> -Cm^l'^ J \W{x)\^/'^dx-Ca^ J \W {x)\'^ dx - C ma-^ 

M 

- C^{m^/^{ma-y/'^ + a^{ma-^)'^)\{x\lam-^ < \x - Rj\ < |am-^}| (166) 

M 

+ ^ TV [ej (V-a-2A + m2a-4 - ma'^ - Cma'^ - _^ - Ci/ma"^) Oj] . 
By the translation invar iance of —A, the last term equals 



MTr[e(V-a"^A + m?a-^ - Cma'^ - Cvma'^ - iy\x\-^)e] _ , 
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and using the Lieb-Yau inequality (and the properties of 6 and that i^a < 2/7r), 



> a-^TY[e{\/^-'^\x\-^ -Cma-^)e]_> -Cma-'^ . (167) 
Further, by the assumption on W , and the assumption min^^^ \ Rk — R(\ > 2am~^, 
j \W{x)\^l'^dx 



am~ 1 / 4<(iR. {x)<oim~ i 

M 



am~^ /4:<\x—Rj\<am~^ 



and, since va < 2/7r, 



J \W{. 

am~ 1 /4<d-[i {x)<am~ ^ 
M 



a;)|^ dx 



< Cv'^c? f {\x - Rj\-^ + a-^mj^dx < Cv'^c?m < Ciy^/^a^^^m . 

■' am-i/4<|3;— -Rj|<am-i 



It follows from this, (167), and (166) that 

TrfV-a^^A + m?a-'^ - ma~'^ + W{x)]- 

> -Cma-'^ -Cv^l'^Q^I'^m-Cm^l'^ j \W{x)\^/'^ dx - Ca^ j \W{x)\'^dx. 

Since ma"^ < Cv^^^^a^^^m when G [3/(167rM), 2/7r] this proves (32) in this case. This 
finishes the proof of Theorem 16. □ 



Proo/ o/ r/ieorem i 7 (Correlation inequality). The proof essentially uses superhar- 
monicity and positive definiteness of |a;|~^ (sec [19]). By superharmonicity and the spherical 
properties of $s. 



- y\-^ > J ^s{z -x)\z- z'\-^ $,(/ - y) 



dzdz' . 
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Also note that for the Coulomb energy, D{^s,^s) = = Cs^^. Therefore, we 

immediately get that 

-1 



E 



I Xj^ ^ j I 



l<i<j<N 



> ^ / — Xj) Iz — z'l ^ ^s{z' — Xj)dzdz' 
l<i<j<N 

$5(2; — Xj)^ \z — z'l'^ ^2 ^s{z' — Xj)^ dzdz 

<3<N l<j<N 

-f y ^s{z)\z- z'\-^^s{z')dzdz' 
= hfiH Mz-Xj)-piz))\z-z'\-'{ Yl Mz'-Xj)-piz'))dzdz' 

i<j<N l<j<N 

+ / p{z) \z - z'\~^ ^ ^siz' - Xj) dzdz' 
i<j<Ar 

-^j p{z)\z - z'\-^ p{z')dzdz' - ND{^„<^s) 

> {p*\x\~^*^s){xj)-D{p)-CNs-'^. 

In the last inequality we have used the positive definiteness of \x\~^ and dropped the first 
term. This proves inequality (33). □ 

Proof of Corollary 21 (Estimate on p^^ * [-"^^l ""^ * (^0 ~ ^t))- Let, with dr as in (4), 

g^{x) = min{(ir(a;)-^/2^ dr(x)-^} . (168) 

We claim that for some constant C > 0, 

|Vp^^* < Cg,:{x). (169) 

To prove this we distinguish two regions. 

First, let dr{x) > 1. Then gr{x) = drix)'"^. Using that |Vy'^^|(x) < Cgrixfdr{x)-'^ = 
drix)"^ we obtain 

M 

M 

< Zj\x — rj|~^ + dr{x)~^ < Mmax{zj}( min \x — rj\) + dr{. 

= Cdr{x)-^ + dr{x)-^ < Cgr{x) . 
If on the other hand dT{x) < 1, then, by using (36) and (41), 

|VpT^*|^r'| < j \Vp"'''{y)\\x-y\-Uy<C j g,{yfdr{y)-^\x-y\-Uy 

M 

= j 9r{.yfd,{y)-^\x-y\-^dy. 



X) ' 



<ir{y)=\y-rj\ 
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Since griy) < dr{y) = \y — fj\ ^^"^ when dr(y) = \y — rj\ this imphes that 



M 



|vpT%N-i| < cf2 f\y-rj\-'^'\^-yr'dy 

M 

= C^\x- rj\-^/^ < CMdrix)-^/^ = Cgr{x) . 

This finishes the proof of (169). 

Let us now proceed to prove inequahty (43). That the difference is positive is again just 
superharmonicity of |a;|~^. It is easy to see that 

\dr{x)-dr{y)\<\x-y\. (170) 

In the case when dr{x) > 2t we can conclude that 

p^^ *\x\~^ ~ p^^ *\x\~^ *(tt < sup { |Vp^^ * } \x - y\ ^t{x - y) dy 

\z-x\<t J 

< Ct sup gr{z) < Ctgr{x) . 

\z-x\<t 

In the last step we have used that if dr{x) > 2t and |z — x| < t (this condition stems from 
the support of $t), then inequality (170) guarantees that ^dr{x) < dr{z) < ^dr{x). This in 
turn implies {^)-'^gr{x) < gr{z) < {ly^f^grix). 

If, on the other hand, dr{x) < 2t, then we claim that 



\p'^^*\x\-'-p'^^*\y\-^<C\x-y\'/\ 
This can be seen as follows. 



(171) 



|pTF*|a;|-l-pTF^|y|-l| 



£f^^,^^.^exHi-o)y\-^). 

^ V(^p^^ *\9x + il~e)y\-^yix-y)d9 

< C f gr{ex+{l-e)y)\x-y\de 

Jo 

M 1 

< C\x-y\y2 \ex + {1 - e)y - rj\-^/^ de 

.7 = 1 "^0 



M ^1 

= Clx-y^Yl 



G{x - y) ^ y - rj 



\x-y\ \x-y\ 



-1/2 



d9. 



Let n = {x-y)/\x-y\, b = {y-rj)/\x-y\, and c = n-b. Then \dn+b\^ > \e+n-b\'^ = {O+cf 
Therefore 



I TF I 1-1 TF 

\p * \x\ " 



M „i 



pTF * \<C\x- y| y^f\e + c\-^l^ dO . 
The integral \0 + c|~i/^ dO is bounded uniformly for c G M. This proves (171). 
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This allows us finally to show that for dr{x) < 2t, 

p^^ * \x\~^ — p^^ * \x\~^ * = J [p^^ *\x\~^ — p^^ *\y\~^)^t{x — y)dy 

< C J \x- y\^/^ -y)dy = Ct^/^ . 
This finishes the proof of the corollary. □ 

Appendix B. Estimates of semi-classical integrals 

In this appendix we give the remaining arguments on the analysis of the integrals in the 
semi-classical proofs of Lemma 33 and Lemma 34. 

Proof of Lemma 33 (Lower bound on Tr[0ff^0]_): Estimate of integral (131). It re- 
mains to estimate the integral in (131). Note that by Taylor's formula for a we have 

^SKp) > a{v,p)+l,,p{u-v,q-p)-C\u-v\{h-^ + \u-v\'') (172) 
-C\q-p\{h-^ + \q-p\^), 

where 

lv,p{'^, q) = is^'^('"^P) - (1 - XI ^A'^i3ip)(iiQj ~ hYl 9idjV{v)uiUj . 

We have used that \Aa{v,p) — Aa{u,q)\ < C\u — v\ + C\q — p\, and similarly, when replacing 
didjF{q) by didjF{p), and didjV{u) by didjV{v). We get that: 

Hi%v,p)<0 => \p\<C{l + \u-vf + \q-p\^). (173) 

(Note that this holds also for e = 0, and uniformly in /3 G [0, 1]; to see the latter, use that 
Tpip) > Ti{p).) This implies that 

j {\u-v\'^ + \q-p\^)Gi,{u-v)Gb{q-p)dpdudq<Ch-'^/'^ , (174) 

and 

j {\u-vr + \q-pr)Gb{u-v)Gbiq-p)dpdvdq<Cb-"'/^ . (175) 

M:i(f,p)<o 
From this we obtain that 

J Gbiu - v)Gbiq - p) \Hi'l{v,p)\ Jpdqdv > - G, (176) 



and hence from (131) that 

Tr\cbH^(b]_ 



> J cP{v + h'abiu-v)fGb{u-v)Gb{q-p)[H^,q{v,p)]_-^^^ 

U&B2 

Gh-'^ib-^'^ + h%) . 



dudq , , 

dvdp 
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Here wc have used the fact that the ^-integration is over a bounded region. From now on 
we may however ignore the restriction on the u-integration. We note that, by using (173) 
and (174), that (f) has support in Bi, and that 6 > 1, we get that 

J (t){v + h^ab{u - v)f{\u - v\{b~'^ + \u-v\^) + \q-p\{b~'^ + \q-p\^) 



hL',1(v,p)<0 



X Gb{u- v)Gb{q - p) dudqdvdp < Cb 



Using this and (172) we find after the simple change of variables u ^ u + v and q ^ q + p 
that 



TV[<^lf^</.]_ 
> J (t){v + h^abufGbiu)Gbiq) 



X [a{v,p) + e.,p(n, q) - C\u\{b-' + \u\'') - C\q\{b-' + \q\')] _ ^^dvdp 

- ch-'^ib-^/^ + h'^b) 

(j){v + h^abufGb{u)Gb{q) \^{v,p) + ^v,p{u, q)] _ dvdp 

- Ch-'\b-^''^ + hH) . (177) 

At this point we divide the (i;,p)-integration into three regions given in terms of a parameter 

A > by 

^- = {{v,p) I ^{v,p) < -A} , n+ = {{v,p) I a{v,p) > A} , J^o = {{v,p) \ \^{v,p)\ < A} . 

The parameter A will be chosen such that 1 > A > Cb~^ for some sufficiently large 
constant C. This is possible if r is small enough and hence b large enough. Then, since all 
the second derivatives of a are bounded we may assume that ^\Aa{v,p)\ < A/2 for all 
{v,p), uniformly in (3. 

We first consider ri_|-. We see from (177) that we only need to integrate over the set 
{iu,q) I C(|itp + |gp) > A}. Also notice that a{v,p) > ^\p\ - C (since T^(p) > Ti{p)) shows 
that we only need to integrate over the set {p \ \p\ < C(l + + |wp)}. Therefore, 

J (t){v + h^abufGb{u)Gb{q) \o{v,p) + iv,p{u, q)] _ dudqdvdp > - Ge~'^^^ . 
{v,p)en+ 

A similar argument shows that on 0_ we can ignore the negative part [ ]_ paying the 
same price — C/i^^e^*^^"^. 

For {v,p) G Q- we estimate the integral 

J (t){v + h^abufGb{u)Gb{q)[a{v,p)+^y,p{u,q)] dudq > {(l){vf + Gh^b)a{v,p) - Ch^b . 

Here we have expanded cfP' to second order at the point v and used the crucial fact that 

e.,p(^i, q)Gb{u)Gb{q) dudq = . (178) 



{v,p)eQ-un+ 
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For {v,p) G r2_ we have, of course, a{v,p) = a(v,p)-. Since the volume of is bounded 
by a constant we get for the integration over ri_ U , 

J (f>{v + h'^abu)^Gi){u)Gb{q) [a{v,p) + ^v,piu, q)] _ dudqdvdp 

> j (f){vfa{v,p)^dvdp-C{hH + e-^^^). (179) 

Finally, let {v,p) G Qq. Observe that, with = {2t + pt'^)'^/^, 

J dp = c„(i?(-[A - Viv)U) - ^(-[A + V{v)U)) < CA , (180) 

by the mean value theorem (uniformly in v). Now, 

(t){v + h^ahuf Gb{u)Gb{q) [a{v,p) + ^„,p('u, q)] _ 
> (t){v + h^abufGb{u)Gb{q)[a{v,p)]_ 

- C(j){v + h'^ahuf Gb{u)Gb{q){b-^ + |np + |gp) , 

and, using the observation above and making the change of variables v ^ v — h^abu in the 
i;-integral, 

J (t){v + h^abuf Gb{u)Gbiq)ib'^ + \uf + \q\'^) dvdpdudq < CAfe"^ . 
{v,p)eQo 

Expanding (fy^ to first order at v we have that 

J (f){v + h'^abuf' Gb{u)Gb{qy5{v,p)- dvdpdudq 
(i;,p)eno 

> j (f){v)^a{v,p)- dvdp + Ch^ab J \u\Gb{u)Gb{q)d-{v,p)- dudqdvdp 
{v,p)efio (v,p)eCio,vesuppV 



> J (j){vfa{v,p)-dvdp-Chb^''^k^. 



As a consequence, 



'{v + h^abu)^ Gb{u)Gb{q) [a{v,p) + (,v,p{u, q)] _ dudqdvdp 

{v,p)eno 



> j<p{vfa{v,p)-dvdp-CK{Khb^/'^ + b-^). (181) 



{v,p)&Vlo 
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Since 

^2 



(v) a{v,p)- dvdp 

= (1 - e)-" J <j){vf [ V + (1 - e?l^-^ - (1 - e)f3-^ + V{v)\ _ dvdp 

> (l-e)-"- J (j){vfap{v,p)-dvdp 

> J (j){v)'^a/3{v,p)^ d.vdp — Ce , 

the lemma follows from (177), (179), and (181) if we choose b = h'^/^. □ 

Proof of Lemma 34 (Construction of a trial density): Estimates of integrals. We 
give here the remaining arguments on the analysis of the integrals in the semi-classical 
proofs of Lemma 34. 

The energy: proof of (132). It remains to estimate the integral in (139). 

Using (173), that Tf^(p) < ^p^, and that hu,q{v,p) = unless u G i?2, we get that 



/ 



x[K,q{v,p)] Gb{u - v)Gb{q - p)(l + Tp{p + h^ab{q - p))) dudqdvdp < C . 
This implies that 

< J x[hu,q{v,p)]Gb{u-v)Gbiq-p)4>{v + h'^ab{u-v))'^ 

X a(v + h^ab{u -v),p + h^ab{q - p)) .f^f^ dvdp + Ch^bh^"' . 

(ztt/i)'* 

From (137) we may now conclude that 

< J xHv,p) + q) - C\u\ib-' + \u\^) - C\q\ib-' + \q\^)] Gb{u)Gb{q) 



UEB2—V 



X (l){v + h^abufaiv + h^abu, p + h^abq) /^f^ dvdp + ChHh-'' . (182) 

{2Trh)'"' 

At this point we introduce the same partition of the (u,p)-integration into sets ft±, J7o as in 
the proof of the lower bound above (with e = 0) with the same A = 6-1/2 ^ /^2/5, 

Then for the integration over ri_|_ we have as above that C(|'up-|-|g|^)>A and hence 

J xHv,p) + C.,p(u, q) - C\u\{b-' + \uf) - C\q\{b-' + \qf)] 

X ^(u + h^abufa{v + h^abu,p + h^abq) Gb{u)Gb{q) dudqdvdp < C e~'*^ < ChH, 

where we have used (136) and that 4> is supported in the ball Bi. 

Similarily, if {v,p) G r2_ then for the (u, g) -integration we can safely assume that the 
argument of x is negative to the effect of paying the same e"*^*^ price. Likewise we may 
ignore the restriction u & B2 — v, since u ^ B2 — v and v + h^abu G Bi implies \u\ > 
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(1 — h^ab)^^ > 1. Expanding 0^ and a to second order at {v,p) G 0_ and using the fact 
that all their second order derivatives are bounded together with (130) we get that 

J xHv,p) + ^v,p{^,q) - C\u\{b-' + \u\^) - C\q\{b-' + \q\^)] 

X (f)(v + h^abu)^a(y + h^abu,p + h^abq) Gi,{u)Gi,{q) dudq 

< J {(j){vf + h^abu ■ V((/)^)(i;)) {a{v,p) + h^ab {u, q) ■ Va{v,p)) Gb{u)Gbiq) dudq 

< (l){vfa{v,p)^+Ch%. 

It is important here that a and Vo" are bounded uniformly in /? < 1 on This follows 
from (136) and jVcj/3(u,p)| < C{1 + \p\). Indeed, (136), in particular, implies that fl- is a 
bounded set (uniformly in (3). The fact that the volume of Q- is bounded also gives that 
the contribution from r2_ to the integral on the right side of (182) is bounded above by 



(27r/i)-" / 4>{vfa(v,p)-dvdp + Gh%h-''. 



Finally, we consider {v,p) G J7o- If we expand to first order at v and a to second order 
at {v,p) and use that all second order derivatives of a are bounded and that Va{v,p) is 
bounded for {v,p) € 0,q we obtain that 

(l){v + h^abufa{v + h^abu,p + h^abq) < (t){vfa{v,p) + Ch^ab{\u\ + \q\) + Ch^aH'^{\u\^ + \q\^). 
This together with the estimate \x{x + y)x — x{x)x\ < \y\ implies that 

/ xHv,p)+^vA^,q) - C\u\ib-^ + \uf) - C\q\{b-^ + \qf)] 

X (l){v + h?abu)'^a{v + h?abu,p + h^abq) Gi,{u)Gh{q) dudq 

< J xHv,p) + ^v,p{^, q) - C\u\ib-' + \u\^) - C\q\{b-' + \q\^)] Gb{u)Gb{q) dudq 

x4>{vfa{v,p) + Ch^b^^^ 

< ^{vfa{v,p)- + C{b-^ + h%^/^) . 

We have here again used that the effect of removing the restriction u ^ B2 — v causes a 
smaller error than the last term above. Note that u £ B2 — v and v + h?abu G Bi imply 
l^^l < 3(1 — h^ab)~^ < 6 and hence we only need to consider \v\ < \v + u\ + \u\ < 8. If we 
use that (180) implies 

Vol(no n {v I \v\ < 8} X W) < CA 

wc see that the contribution from to the integral on the right side of (182) is bounded 
above by 

(27r/i)-'^ / (P{vfa{v,p)- dvdp + Ch'^'ib-^ + h'^b^/^)A . 
This finishes the proof of the upper bound on the energy in (132) . 

The density: proof of (133) and (134). Here it remains to estimate the integral in 
(140). The strategy is to freeze the variable \p\ in ^y^p so that the remaining dependence 
on \p\ is explicitly integrable. This is accomplished in the estimate (183) below. After the 
|p|-integration the proof is almost the same as in the non-relativistic case [36]. 
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We write p = \p\(jj and define 

Po = {f3\V{v)-f + 2\V{v)-\)'/^u = r7(F(i;))V"a; . 
We will then prove that if u e B2 — v then 

x[a(v,p)+^v,poiu,q) + R{u,q)] < x[hu+v,q+p(.v,p)] < xW{v,p) + ^v,po{u,q)-Riu,q)] , (183) 
where 

R{u, q) = C{\u\{b-^ + \uf) + i\q\ + A)(6-^ + + (^^ + \uf + + |g|')A-') . 

We first prove (183) for {v,p) G Qq. In this case wc have 

r]{V{v) + A)2/" < _p2 < ^(y(^) _ A)2/» , 

from which it follows Ihrit — Pq\ < CA with a constant independent of /3 G [0, 1]. 

Let G{t) = ^/iFH^W^ - [i-^ , so that Tgip) = G{p^) (we suppress that G depends on 
(3). Note that then didjTpip) = ApipjG" (p"^) + 26ijG'{p'^), and so, in particular, AT(j{p) = 
4G"{p^)p'^ + 2nG'{p'^). Therefore, using that pi = \p\u}i,po^i = \po\ijOi, 

\iv,p{u,q) - ^y^PQ{u,q)\ 

< ^ \Aa{v,p) - Aa{v,po)\ + ^ ^ \didj[Tf,{p) - Tp{po)]\ Wiqj] 

< C{b'' + \q\'){\G"{p')p' - G"{pI)pI\ + |G'(/) - G'ipDl) 

< C{b-^ + \q\'^)\p'^ -pI\ < CA(6-V|g|^). (184) 

Here we have used the choice of pQ and that G'{t) and tG"{t) have bounded derivatives 
uniformly in /? € [0, 1]. If we combine (137) with (184) wc obtain that 

\K+v,q+p{v,p) - aiv,p) - ^v,poiu,Q)\ < C\u\{b-^ + |u|2) + \q\{b-^ + \qf) + CA{b-^ + \qf) , 

which is, in fact, stronger than (183). 

If {v,p) £ we see that the left inequality in (183) is only violated if Cv,poiu,q) < —A 
and the right inequality is only violated if hu+v,q+p{v,p) < 0. Since {v,p) G we must in 
both cases have C(|up + Iq]"^) > A. We hence get (again using (137)) that 

\hu +v,q+p 

< \hu+v,q+p{v,p) -ct{v,p) -Cu,p{^^Q)\ + \^v,p{^^l)\ + \Cv,pQ{u,q)\ 

< C\u\{b-^ + |u|2) + G\q\{b-^ + \q\^) + C{b-^ + \u\^ + |g|2) 

< C\u\{b-^ + \uf ) + C\q\{b-^ + |g|2) + C{b-^ + + \q\^){\u\^ + \q\^)A-^ , 

which gives (183) in this case. 

Finally, if {v,p) £ Q- then the left inequality in (183) is only violated if hu+v,q+p{v,p) > 
and the right inequality is only violated if ^v,po{u,q) > A. In both cases this implies that 
C(|tip + |gp) > A and hence the same argument as for proves (183). 

Using (183) we can estimate the density in (140) from above and below. We will discuss 
the lower bound on the density. The upper bound is similar. Performing the |p|-integral in 
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(140) we obtain 

p^{x) > J E{u, q, V, u)Gb{u)Gb{q)G(^h%)-i {x - v - h?abu) dvdw (185) 



u^B2—v 

E{u, q,v- h'^abu, uj)Gb{u)Gb{q)GQ^2b)-i {x - v) dvduj ^^^j^y , 



I 



where 

3{u,q,v,Lo) = n'^ri {V {v) + i^^p^{u,q) + R{u,q)) . 

Wc have that 



V{v - h abu) + iv-h^ahu,po{'^, q) + R{u. q) 
< V{v) - h^abuVV{v) + Cv,po{u, q) + R{u, q) + Ch'^a^b^\u\^ + Ch^ab\u\{b-^ + \u\ 



2\ 



\ri{s)-ri{t)-ri'{t){s-t)\ < 



< V{v)-h''abuW{v)+S,^^p,,{u,q)+R{u,q)+Ch'^a'^b'^\u\\ (186) 

(In the last hne we have used that h^ab < 1 and the definition of R{u, q).) We now use that 

C|s-t|3/2 + C(|s| + |t|)|s-t|2, n = 3 

C(|s|5-2 + |i|f-2 + |s|-2 + |t|«-2)|s_i|2^ n>4 • 

(187) 

We continue with the case n = 3 and leave n > 4 to the reader. 

If we use the fact that r]'{V{v)) is bounded independently of /3 G [0, 1] we obtain from 
(186) and (187) 

nE{u,q,v - h^abu,uj) > r]{V{v)) + r]'{V{v)){^v.p^,{u,q) - h'^abu\/V{v)) 

- C{b-^ + \q\^ + \u\^ + h^ab\u\ + R{u, q)f'^ 

- C{b-^ + \qf + |u|2 + h^ab\u\ + R{u, q)f 

- CR{u, q) - Ch'^a%'^\u\^ - Ch^ab\u\ . (188) 
It is now crucial that (see (130) and (178)) 



/ 



{^v,po{u, q) - h^abuW{v)) Gb{u)Gb{q) dudq = . 

,2 



Since v G supp(y) C B2/2 and (1 — h ab)u ^ B2 — v imphes \u\ > 1/2 we find 



/ 



<Ce-^/^ <Ch^/^ . (189) 



i^v,poiu, q) - h^abvS/V{v)) Gb{u)Gb{q) dudq 

{l-h?ah)ueB2-v 

Combining this with (188) and inserting it into (185) we arrive at (recall that A = b~^/'^) 

p^{x) > (27r/i)"^W3 j ri\y{v)\GQ,2b)-i{x - v)dv 

- Ch-^ (/i^/^ + 6-3/2 ^ (/,2^^^3/2^-3/4) _ (^qq^ 

Here we have removed the constraint (1 — h'^ab)u ^ B2 — v hy the same argument as above. 

We shift the i;-coordinate by x, and then expand r7[F(x + v)^ in the integral at x by 
expanding V to second order at x and using (187). Then 

r] [V{x + v)] >r] [V{x)] + 77' [V{x)]VV{x) ■ v - C{\v\^/^ + \v\^)- 
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Then we obtain from (190) (using (130)) that 

p^{x) - {2Trh)-'^uJs7][V{x)] > -Ch-^{h^/^ + (hHf*) > - . 

This finishes the proof of (133). 

Lastly, we prove (134). By integrating (190) we see that 

j 4>{x)'^ p^{x) dx 

[V{v)] dxdv - Ch-^^^/^ 

> (27r/i)-3w3 j 4>{vf T] [V{v)] dv - Ch-^+^/^ . 
In the last step we have expanded to second order at v to obtain that (see also (130)) 

j (t>ixfG(h'^t,)-^ix -v)dx< cp{vf + Ch^l^ . 

This finishes the proof of (134) and therefore of Lemma 34. □ 
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